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Abstract: We analyze transverse thrust in the framework of Soft Collinear Effective 
Theory and obtain a factorized expression for the cross section that permits resummation 
of terms enhanced in the dijet limit to arbitrary accuracy. The factorization theorem for 
this hadron-collider event-shape variable involves collinear emissions at different virtualities 
and suffers from a collinear anomaly. We compute all its ingredients at the one-loop order, 
and show that the two-loop input for next-to-next-to-leading logarithmic accuracy can be 
extracted numerically, from existing fixed-order codes. 
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1 Introduction 

Event-shape variables are designed to measure geometrical properties of energy flow in 
collider events. They were among the first observables proposed to test quantum chromo- 
dynamics (QCD), and can also be used to discriminate beyond-the-Standard-Model (BSM) 
physics against the QCD background. Numerous event-shape studies have appeared over 
the years, notably including extractions of the strong coupling, a s . The majority of the 
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existing work has focused on leptonic collisions, and on deep-inelastic scattering (DIS). 
Nevertheless, event shapes are also of great interest in the much richer environment of 
hadronic collisions. Indeed, a lot of recent work is using event shapes as a tool to study 
jet substructure. Event shapes could also be instrumental to get a better handle on some 
poorly understood aspects of hadronic collisions, such as the underlying event. 

Several event-shape variables for hadronic collisions were studied in refs. [1, 2], They 
were defined in analogy to the ones in leptonic collisions, but in terms of the components 
of the three-momenta transverse to the beam direction. In this paper we focus on the 
archetypal event shape, thrust. Transverse thrust, which we denote by Tj_, is defined, in 
analogy to ordinary thrust, as 


Tj_ = max 

n± 


I Pi± • 

£t Iml 


(i.i) 


where the sum is over all the particles in the final state, with momenta p). Throughout the 
paper, the subindex T denotes the (two) momentum components transverse to the beam 
direction. The vector n± which maximizes the ratio on the right-hand side of eq. (1.1) is 
called the transverse-thrust axis. Transverse thrust has been measured at the Large Hadron 
Collider (LHC) [3-6] and previously also at the Tevatron [7]. Here, we will study this 
quantity in the dijet limit, where Tj_ —> 1, and obtain a factorization formula that allows 
us to resum the enhanced terms arising in this limit. Resummation for hadron-collider 
event shapes in the dijet limit at next-to-leading-logarithmic (NLL) accuracy was studied 
in refs. [1, 2] within the automated resummation framework CAESAR [8]. In the present 
paper we analyze transverse thrust using Soft Collinear Effective Theory (SCET) [9-11] 
(see [12] for an introduction), and obtain an all-order factorization formula that allows for 
resummation at any desired accuracy. 

To derive the factorization theorem, we start with the lepton-collider case. There is 
no need to restrict oneself to the plane transverse to the beam in this case, but doing so 
provides us with a simpler environment to analyze the factorization of transverse observ¬ 
ables. In contrast to standard thrust, Tj_ —> 1 does not imply that the event consists of 
two low-mass jets. Nevertheless, the terms which are enhanced in this limit do arise from 
two-jet configurations, such as the one shown on the left-hand side of figure 1, in which 
all the radiation is soft or collinear to the two low-mass jets. The resulting factorization 
formula has the same structure as the one for thrust. It involves a hard function which col¬ 
lects the virtual corrections to the hard scattering process, two jet functions describing the 
collinear emissions and a soft function. The energy of the soft emissions is parametrically 
lower than the typical mass of the jets. 

In the hadron-collider case also the incoming partons carry color charge and the ef¬ 
fective theory involves additional collinear fields which describe the initial-state radiation. 
The proton matrix element of these fields defines beam functions, which can be factorized 
into a perturbative kernel, describing the emissions, convolved with the standard parton 
distribution functions (PDFs). Interestingly, the virtuality of the initial-state collinear 
fields is parametrically of the same order as the one of the soft fields, and is lower than 
the virtuality of the collinear fields of the final-state jets. As is typical for problems which 
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Figure 1. Schematic representation of a dijet event in the T± —> 1 limit for leptonic (left panel) 
and hadronic (right panel) collisions. The soft radiation s and the collinear emissions ci, C 2 , c a , Cb 
are represented by different fields in the effective theory. The typical virtuality of the fields c a , Cb 
and s is the same, and is lower than the virtuality of c\ and C 2 . 

involve soft and collinear fields of the same virtuality, transverse thrust suffers from a 
collinear anomaly: the soft and beam functions are not well defined individually and their 
product involves large logarithms associated with the large rapidity difference between the 
emissions from the two incoming particles [13]. To compute the beam and soft functions 
individually, one needs to introduce an additional regulator, which can be removed after 
combining the functions. Traditionally, this regularization was achieved by taking the Wil¬ 
son lines describing soft and collinear emissions in these functions off the light-cone, see 
e.g. [14], However, in an effective theory context, it is more convenient to use an analytic 
regulator which does not introduce additional scales into the problem. The cancellation of 
the divergences in the additional regulator imposes constraints on the form of the large log¬ 
arithms generated by the collinear anomaly. These constraints are particularly interesting 
in our case due to the nontrivial color structure and angular dependence of the soft func¬ 
tion for transverse thrust. The fact that the problem involves nontrivial color structure, 
collinear fields at different virtualities and a collinear anomaly illustrates that factorization 
for transverse thrust is quite nontrivial. 

The resummation of large logarithms is achieved by solving the renormalization group 
(RG) equations of the ingredients in Laplace space. Transforming back to momentum 
space, we provide an analytic form of the resummed partonic cross section. Towards 
the goal of achieving next-to-next-to-leading logarithmic (N 2 LL) accuracy we evaluate all 
the constituents of the theorem at one-loop accuracy. The other ingredients for N 2 LL 
resummation are the two-loop anomalous dimensions and the two-loop anomaly coefficient. 
Using factorization constraints, we show that the only unknown quantities are three two- 
loop coefficients. We determine one of these coefficients numerically by comparing to the 
next-to-next-to-leading order (N 2 LO) fixed-order result for transverse thrust in leptonic 
collisions. This also provides a numerical check on our factorization formula. We then 
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show that also the remaining coefficients can be extracted numerically, by considering 
transverse thrust in Drell-Yan and Higgs production events, and using existing N 2 LO fixed- 
order codes for Higgs and Drell-Yan production. This determination, together with a 
numerical implementation and phenomenological analysis of the resummed cross section, 
will performed in a future publication. 

Our paper is organized as follows: We derive the factorization formula for the transverse- 
thrust differential distribution in section 2 and compute its ingredients at the one-loop level 
in section 3. Solving the associated RG equations, we derive in section 4 resummed expres¬ 
sions for the cross section in the dijet limit. In section 5, we first compare our resummed 
results for the lepton-collider case with the fixed-order computation. This allows us to de¬ 
termine the two-loop anomalous dimensions, so that we have all the ingredients for N 2 LL 
resummation in the lepton-collider case. We then show that a similar procedure can be 
used in the hadron-collider case. Our conclusions and an outlook on future work are pre¬ 
sented in section 6. The appendices collect anomalous dimensions, and provide details on 
the one-loop computations of the jet, soft and beam functions. 

2 Factorization formula 

Our goal is to derive a factorization formula for the transverse-thrust differential distribu¬ 
tion that is valid in the dijet limit. If we define, as usual, 

tj_ := 1 — Tj_, (2.1) 

the dijet limit, where the event contains two low-mass jets, corresponds to t± —> 0. However 
unlike the usual thrust case, the limit r± —> 0 contains not only dijet configurations but also 
configurations where all the particles lie in a plane which contains the beam, see figure 2. 
The dijet configurations give singular perturbative contributions to the cross section at low 
t± 

da a s 
dr± Tj_ 

while the multi-jet configurations are regular and thus power-suppressed at low tj_. In 
the following, we will study the singular terms in the limit t_l —> 0 in detail, and resum 
their contribution to all orders in perturbation theory. The power-suppressed terms can be 
added by matching to fixed-order results. We will find that the matching corrections are 
larger for transverse than for regular thrust, the reason could be the presence of multi-jet 
configurations even at low tj_. 

In order to understand how one should treat event-shape variables that involve only 
momenta transverse to the beam direction in SCET, we will first consider leptonic collisions. 
At lepton colliders, there is not much experimental motivation to restrict event shapes to 
the transverse plane, but doing so provides us with a simplified environment to analyze 
factorization for such observables. A typical dijet configuration in a leptonic collision is 
depicted in the left panel of figure 1. The incoming electron and positron are taken in the 
^ direction, and are represented by the black arrows in the figure. The jets are emitted 
at an angle 8 with respect to the beam. As for the usual lepton-collider event shapes, 
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Figure 2. Left: A two-jet configuration with low r_j_. The figure shows the thrust axis n (green), 
the transverse thrust axis n± (red), and the beam as a dashed line in the z-direction (blue). Right: 
A planar three-jet configuration with t± = 0. Any distribution of particles that is restricted to a 
plane which contains the beam has t±_ = 0. 

the effective theory needs two collinear modes, which we call ci and C 2 , to account for the 
energetic particles in the final state (represented by the red lines in figure 1), and a soft 
mode s which describes final-state soft radiation (represented by the orange wavy lines in 
figure 1). We will perform the detailed factorization analysis below, but the astute reader 
will have guessed that the result will be a factorization formula of the form 

« H{Q) ■ (Jc 2 ± ® Jc i_L <8> S±)(t±) . (2.3) 

The hard function H(Q) collects the virtual corrections to the hard-scattering process and 
is the same as for regular thrust. The jet and soft functions, on the other hand, differ from 
the standard case, because the phase-space constraints associated with t± only act in the 
transverse plane. 

The lepton-collider case can be contrasted with the hadronic case depicted in the right 
panel of figure 1. A complication on the kinematic level is that the partonic collisions, 
which produce the jets at hadron colliders, are not taking place in the hadronic center-of- 
rnass frame and do not have fixed energy. Therefore the jets are not back-to-back and their 
energies are not fixed. To be able treat the process perturbatively, one needs to ensure 
that each event involves an underlying hard collision. One way to do this is to impose a 
minimum transverse momentum. 

A complication for the theoretical description is that there is also initial-state radiation 
and the effective theory includes two additional collinear modes c a and q, along the beam 
directions. The initial-state radiation is described by beam functions, which are proton 
matrix elements of these collinear fields. For perturbative values of tj_, these can be 
factorized into PDFs convolved with perturbative kernels describing the emissions. We 
will show below that the collinear modes c a and Cb have a low virtuality, of the same 
size as the one of the soft fields, and that the convolution of the beam functions and the 
soft function suffers from a collinear anomaly [13], i.e. the functions are not individually 
well-defined and their product has a logarithmic dependence on the hard scale. 
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2.1 Leptonic collisions 


Let us now analyze the leptonic case in detail. In the dijet limit the final state X consists 
of particles X C1 and X C2 with large energies flying along the jet direction, defined by the 
thrust axis n. and soft particles X s which are radiated at arbitrary angles. Up to power 
suppressed terms, we can thus simplify the denominator in the transverse-thrust expression 
as 


r, \Pi±\= Ei \ sin#j| ~ | sin0| ^ Ej + ^ | sinO^Ei ~ Q|sinfl|, (2.4) 

*Sci,C2,s ieci,C2,s *Gci,C 2 iGs 


where Ej is the energy of the particle with 3-momentum p), and Q is the total center-of- 
rnass energy. The sum over i in eq. (2.4) runs over the three sectors in the effective theory. 
In the second step we used the approximation that the collinear particles fly approximately 
along the jet direction 6j ~ 9 , and we have neglected the small contribution of the soft 
particles in the final step. In order for this approximation to be valid, the angle 9 cannot 
be too small (the scaling of the fields derived below will yield the condition 9 > t±, which 
is always satisfied in the dijet tj_ —> 0 limit). 

For a final state X in the dijet limit, we thus have 


T±(X) 


1 

Q| sin#| 


^2 dm I ~ |™l - ml) = r±(x ci ) + t±(x C2 ) + t±(x s ) . 
iex 


(2.5) 


Eq. (2.5) is in a suitable form to derive a factorization theorem in SCET, because the 
sum over particles naturally separates into the different sectors of the effective theory, as 
indicated by the right-hand side of eq. (2.5). From eq. (2.5) one can also read off the 
relevant scaling of the different modes. To do so, let us first introduce the light-like vectors 
rij := (1, Hi) and n t := (1, —Hi), where Hi is the direction of jet i. Therefore, in the lepton- 
collider case we are considering in this section, we have n\ = n 2 , hi = 77 , 2 , and the vector 
H\ =: H is the thrust axis. Throughout the paper, we denote momentum components 
transverse to the thrust vector Hi with the subindex T. Table 1 summarizes the notations 
we employ for the different axis and relevant directions. To analyze the scaling, we split a 
generic momentum p^ into components along the jet and a remainder which is transverse 
to the jet: 


11 -11 LL -11 

X n i / \ n i U n i U 

= (m -p)2~ + (n 1 -p)2~+p9 r =:p + -X +p^2~ +p^ . 


( 2 . 6 ) 


Please note that we use the notation pj and pj to indicate quantities transverse to the 
thrust axis. The notation is instead used to indicate quantities transverse to the beam 
axis and we use the notation p 1 ^j for quantities which are transverse to both axes. The 
contribution of a collinear particle to regular thrust is driven by the small light-cone com¬ 
ponent of its momentum, i.e. by p_ = pj/p+ for the ci particles. For transverse thrust, 
we will see that the relevant quantity is pjj/p + which scales in the same way. Eq. (2.5) 
then tells us that the components of the soft modes scale like the small components of the 
collinear modes. The components (p+,p_,pV) of the different momenta therefore scale as 


ci : ( 1 ,t±,i/tl)Q, c 2 : (t±, 1,-/Tl)Q> s : (t ± ,t ± ,t±)Q , (2.7) 
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direction 


vector 


perp. dir. 


beams 

a 

11 

= (1,0,0, 1) 

T 

b 

n b 

= (1, 0, 0,-1) 

jets 

1 

2 

11 

u, 

n 2 

= (1, sin0i, 0, cos$i) 

= (1, — sin 02 , 0, — cos # 2 ) 

T 

thrust 


n 

= (sin 6. 0,cos 6) 

T 

transverse thrust 


rax 

= ( 1,0, 0) 



Table 1. Summary of the notation used in the text for the relevant axes and directions. The 
last column shows the symbol we use to denote the momentum components perpendicular to the 
respective directions. We use the notation pJJj for the component of which is transverse to 
both axes. At e + e~ colliders one has 61=62 = 6 . At hadron colliders, the same is true in the 
center-of-mass frame of the underlying hard collision. 


which is the same scaling that is relevant for ordinary thrust. The associated effective 
theory is usually called SCETp Given this scaling, one can further expand the contribution 
of the collinear particles to transverse thrust by using the fact that \pijj\ <C \n± • Pi±\, and 
write 


r±(X Ch2 )Q _ L = (IPiJ-l - I Pi ■ nil) = Y 

ieci |2 ieci, 2 

~ v 1 ~ V 1 Ms: ~ 1 v /,)217 f2.8) 

.“ 2 \pi ■ nj_ | 2 |^_l| | sin 6>| 2Ei ’ 

ieci,2 *6 ci,2 *eci j2 

where := Q| sin0|. 

To put forward a factorization formula in SCET, we start from the expression for the 
QCD cross section differential in t± 

= ^iY\ M ( e+e ~ - >X )| 2 ( 27r (q - Px)5{t ± - T ± (X)), (2.9) 

where rj_(X) is given by eq. (2.5). At leading order in the electroweak couplings, the matrix 
element squared can be written as 

\M(e + e~ -4 X)| 2 = Y (2.10) 

i=V,A 



where 

3tM = E'S/W f !‘«/( I ). ( 2 - n ) 

are the vector (V) and axial (A) currents, with Yy = 7 ^, and T^ = 7 M q 5 ; a is a color 
index, and / denotes flavor. In the following we will leave the sum over color and flavor 
implicit. The photonic contribution to the lepton tensor L l ^ v is given by 


L 


v 

fils 


2 Q 2 



r)Pl/^P2^ P2/iPlv 

o 5 



L ,l, = 0, 


( 2 . 12 ) 
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with pi and P 2 the electron and positron momenta, respectively, e the charge of the electron, 
and Qf the electric charge of fermion /. 

The steps needed to obtain a factorized form of the cross section in eq. (2.9) above 
are analogous to other SCET derivations that exist in the literature, the main difference 
being only that we want to leave the angle 9 unintegrated. 1 For this reason, we will just go 
through the main points of the derivation here; we refer to ref. [16] for further details. The 
first step is to match the currents onto SCET operators; then we perform held redefinitions 
on the SCET fields to decouple the soft gluons. After the decoupling, the different sectors 
no longer interact and we can write the final state as 

|X> = \X S ) \X C1 ) \X C2 ) (2.13) 

and the transverse thrust constraint in the form 


<5(ti - T±(X)) = S(t±- t±(X s ) - T±(X C1 ) - T l{X C2 )) 


j dT Cl ±d,T C2 ±dT s ± S(t± t s _ i_ t Ci _l t"c2-l) 


X d(r Cl _L - tj_(X c1 ))6(t C2 _l - t±(X C2 ))6(t s ± - t±(X s )) . 


(2.14) 


We then obtain the cross section as a convolution of matrix elements in the different sectors 
of the effective theory 

^ L ^u\C{Q 2 )\ 2 [dT Cl xdT C2 ±dT s ±6(r±-T s ±-T Cl x-T C2 x) 


dr± 2 Q 2 


i=V,A 


X(E(°I Xia(°) \ X ^) <X C2 | xO'(°) 1°) ~ ^( A «)) 

A c 2 

E <°l 4,^(0) I**) (Xc^xiA 0 ) 1°) - r±(X Cl )) 


, A',. 


. w 


x E 0 y c U 0 ) y C i(°) 




Yl(0)Y C2 (0)] kk/ |o)j(r s j_ - r ± (X0) 


X T i,ap T i,a'P' (2tt ) 4 d^(q - PX C1 - PX C2 ~ PX J , 


(2.15) 


where C is the Fourier transform of the matching coefficient from the QCD current to 
the SCET operators. We use latin indices to denote color and greek indices for the Dirac 
structure. The y’s are quark-jet fields in SCET, which include the quark coming out of the 
hard collision and its interactions with collinear particles. The Y’s are soft Wilson lines, 
which encode soft interactions at leading power. As a final step, we want to expand away 
small components in the momentum conservation 5-function and make the dependence on 
the angle 9 explicit. Up to power suppressed terms, we have 


S {li \q-px ci ~Px C2 -px J = 2<5(ni -p Xci ~ Q)S{n 2 ■ Px C2 ~Q)d {2) {Px cl + Px C2 ) ■ ( 2 -!6) 

1 For regular e + e“ event-shape variables, the dependence on 9 has been considered in ref. [15]. 
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To make the dependence on the angle 9 explicit, we now explicitly distinguish the reference 
vector ft in SCET from the thrust axis Ht, which is derived from the particles in a given 
event. We then introduce 

1= (d 3 H 5^ 3 \H — Ht) , (2-17) 

into eq. (2.15). In the effective theory, the thrust axis is given by Ht = Px C] /\'Px <:] | up to 
power corrections. Inserting this into the above equation and using the fact that momentum 
conservation fixes \px ci \ = Q/ 2, we can rewrite it in the form 

Jd 3 H6^ 3 \H-H T ) = {2tt) J d cos 9 (^j 5 (2) (p£). (2.18) 

After these manipulations, the momentum conservation 5 functions only act on a single 
sector and the cross section factorizes into separate collinear and soft matrix elements. The 
collinear matrix elements define jet functions and can be written in the form 

JcMt C i ±) ■= £ <0| x£( 0) \X C1 ) (X cl \xi, a (0) |0) 

X C1 

X <5(t Ci _l - t±(X c1 ))6(Q - hi ■ px ci )d {2) (px c1 t) , (2.19) 

Jc 2 ±(t C2± ) := £ (0| Xc2,a(0) I X C2 ) (X C2 \ X k C2 ,p(0) |0) 

x S(t C2± - t±(X C2 ))5(Q - n C2 ■ Px C2 )d (2) (px C2 t) , (2.20) 



2(2vr) 


ii 

2 


where the collinear-sector transverse-thrust constraint has been expanded according to 
eq. (2.8). The soft matrix element has the form 


S±(r s ±) : = 




r±(X s )). 


( 2 . 21 ) 

Expanding away power suppressed terms in the ^-functions in eq. (2.15), inserting the 
definitions of the jet and soft functions, and contracting the Dirac structure, we obtain the 
desired factorized expression for the cross section 

da 7r N c Q 2 f a 2 n n f 

= 2Q? ^ + C “ > J 


C)(t_l T s _l 'Tc 2 -1- ^ciJ_) ("5:2^)^ci_L (Tci_L)“S'-L(T sJ .), (2.22) 


with H(Q 2 ) := |C(Q 2 )| 2 , and where a = e 2 /(47r) is the fine structure constant and N c the 
number of colors. To include also the contributions from Z exchange, one should substitute 
Q 2 e 4 /(2Q 2 ) [L v + L a ), where L V,A are defined in appendix C. 

We already stressed above that the jet functions that appear in eq. (2.22) are not the 
inclusive jet functions that one needs, for instance, in ordinary thrust. At lowest order, we 
have 

S±{t s± ) = 6(t s ±) ; J c ±( T c±) = <5(t c _l) ; H(Q 2 ) = 1. (2.23) 

We compute all these functions at one loop in section 3. 
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2.2 Hadronic collisions 


Having obtained the factorization formula in the lepton-collider case in the previous section, 
we now move to hadronic collisions, which is the actual situation of interest. The right 
panel of figure 1 depicts a typical dijet configuration in a hadronic collision. The final-state 
jets do not need to be back-to-back in the lab frame, since the momentum fractions of the 
partons entering the hard interaction can be very different. To define the jet axes, one can 
use a jet algorithm. Since our treatment concerns events with two narrow energetic jets, 
any choice of the algorithm will lead to the same jet directions in the limit t± —> 0, up to 
terms which are power-suppressed in this limit. The jet algorithm provides us with the two 
angles 9\ and 0 2 of the two energetic jets, as indicated in Figure 1, and with their energies 
Ej 1 and Ej 2 . Since the jets are massless, these angles are in one-to-one correspondence to 
the rapidities of the two jets. To ensure that we indeed deal with hard collisions, one needs 
to require that the two jets are hard. This can, for example, be achieved by imposing that 

Qj_ := | sin 0\\Ej 1 + \sm6 2 \Ej 2 > Qo , (2.24) 

note that the definition of Q_l now involves 9 i )2 and Ej 12 , as adequate for hadronic col¬ 
lisions. The scale Qo must be large enough so that the soft scale is still perturbative, 
Q_i_t_ l 3> Aqcd- In the limit tj_ —> 0, the two contributions to Qj_ correspond to the 
transverse momenta of the two, approximately massless jets and momentum conservation 
in the transverse plane requires that the two contributions to Qj_ must be equal. Below, we 
will perform a boost along the beam axis to the frame in which the jets are back-to-back 
and their energies are equal. In the limit tj_ —> 0 this is simply the center-of-mass frame 
of the hard scattering after initial-state radiation. 

The modes that we need in the effective theory include the collinear modes ci j2 for 
the final-state jets and the soft mode s that were already present in the lepton-collider 
case. On top of this, we need two additional collinear modes in the beam directions, which 
we denote by c a fi and describe the initial-state radiation. They are represented by the 
green lines in figure 1. To derive the desired factorization theorem, we start by writing the 
expression for the QCD cross section differential in tj_, in a proton-proton collision, 

— E \ M (PP X )| 2 * (2vt) 4 5 (4) (Pa + P b - Px) S(TJ_ - T ± (X))0(Q ± - Qo), 
arj_ x 

(2.25) 

where P a and Pb are the momenta of the protons, px is the total final-state momentum, 
and Ecm is the hadronic center-of-mass energy. 

Before moving to the factorization analysis, it is useful to set up the kinematics. The 
momenta of the protons are given by 

n A 4 

PH = Ecm^ ; P£ = E CM f, (2.26) 

and we assume that the partons which produce the two jets carry fractions x a and Xb of 
the proton momenta. 2 Following ref. [17], we write the total final-state collinear momenta 

2 The quantities x a and Xb are the momentum fractions after initial-state radiation, they are not equal 

to the momentum fractions inside the PDFs. 
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Figure 3. Virtualities of the different modes present in the hadron-collider case. 


in the beam directions as 

Px Ca = (1 — x a)EcM~^~ + b ^_|_ + b a - — , (2.27) 

11 -LI 

71 77 

p» Xcb = (1 - x h )E CM -f + +b b ^f. (2.28) 

The hrst term is the proton remnant, the remainder arises because the leading parton 
radiates into the final state. The momenta of the partons that enter the hard interaction 
are 

p£ = x a E C MY ~ b» ± - b a . (2.29) 

LL -11 

77 r 71 r 

p' b = x b E CM f - K ± (2.30) 

The total final-state momentum px is given by 

Px = Px s + Px ci + Px C2 + Px Ca + Px Cb , (2.31) 

and momentum conservation P a + P b = px then implies the partonic relation 

Pa+Pb= Px s + Px cl + Px C2 • (2.32) 

As in the lepton-collider case, we can simplify the denominator in the expression for 
transverse thrust by dropping power-suppressed contributions 

y \Pix\ = y Ei\ sin 0i | ~ I sin^il y £) + | sin(9 2 | y E { + y E { sin0*, 

i i i=c\ i=C 2 i=s,c a ,Ci) 

~ | sin 6\ | Ej x + | sin 0 2 1 Ej 2 = Q _ L . (2.33) 


We see that the denominator reduces to Q± in the dijet limit. To obtain this result, 
we have used that in the eg 2 sectors the angles between each particle and the beam are 
approximately equal to the jet direction. The contribution from the soft sector is negligible, 
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Figure 4. Feynman diagram for the Born level qq —* q'q' channel. 


as in the leptonic case, and also the c a ,b sectors do not contribute to eq. (2.33) at leading 
power, since sin0j ~ 0 for particles collinear to the beam. As in the leptonic case, we thus 
have 

Q±r±{X) = ^(|pLl| - |n ± -p i± \) , (2.34) 

i 

and the particle sum separates into sums in the different sectors of the theory 

rj_ (X) = r_i_ (X B ) + t_l (X C1 ) + rj_(X C2 ) + r ± (X Ca ) + rj_ (X Cb ). (2.35) 


From this, we can read off the relevant scaling of the different modes. Obviously, the c\ 2 
and s modes have the same scaling as in the lepton-collider case. In order to contribute, the 
transverse components b 1 ' b± of the c a j, modes need to scale like tj_, and therefore the small 
components of these modes scale like . We are thus in a situation where the collinear 
modes c a ^ and the soft mode s have the same virtuality, which is usually called a SCETn 
problem. The virtualities of all the different modes are summarised in figure 3. Since the 
collinear fields in the jets have a large virtuality, the study of transverse thrust in hadronic 
collisions thus involves, quite uniquely, SCETi and SCETn together in the same problem. 
Since we have collinear and soft modes with virtualities that are parametrically of the 
same order, we should expect a collinear factorization anomaly in the effective theory [13]. 
Below, we will find that such an anomaly is indeed present. 

There are many different partonic channels that contribute to the cross section since 
all four partons involved in the hard collision can be quarks or gluons. The different hard- 
scattering channels correspond to different leading-power operators in the effective theory. 
These operators are built from collinear quark fields Xi and collinear gluon fields Aj , and 
involve one field for each direction. The corresponding operators and their one-loop Wilson 
coefficients were given in [18], and recently these results were extended to two-loop order 
in [19]. To keep the following discussion simple, we focus on the hard process qq —» q'q' 
where the outgoing quarks have a different flavor than the incoming ones. For this case a 
basis of the relevant SCET operators is given by 

Oi(x; s, t, u, v ) = Xb(x + tn b )lA A Xa(x + sn a ) x 2 (z + vn 2 )^t A xi(x + uni), 

0 2 (x; s, t, u, v) = Xb(x + tnb)'-ffj,Xa(x + sn a ) X 2 (x + uh 2 )7 / "xi(^ + uhi). (2.36) 


As usual, the operators are smeared over the light-cone direction conjugate to the large 
momentum flow and their contribution to the effective Lagrangian is obtained after con¬ 
volution with the Wilson coefficients: 


AEscet 



Jds dt dudv Cj(s,t,u,v) Oi(xj s,t,u,v ), 


(2.37) 
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where the sum over the different operators I is implied. The Fourier transforms of the Wil¬ 
son coefficients Cj(n a -p a , n b -p b , n\ - pi, n -2 -pi) depend on the large momentum components 
of the collinear fields and are directly related to the relevant scattering amplitudes [20] ; for 
brevity we will not write the arguments explicitly in the following. In our channel only the 
single diagram shown in Fig. 4 contributes to the coefficients at leading order. 

Due to its color structure, it only generates a nonzero Wilson coefficient C \, but at 
higher orders also the second operator will be present. Under renormalization the two 
operators mix, so that the RG equation becomes matrix-valued. As usual, one can perform 
the necessary field redefinitions in SCET to achieve the decoupling of the soft gluons and 
finds that the resulting soft function is a matrix in color space. After the decoupling, the 
final state is decomposed as 


|X> = \X a ) \X Cl ) \X C2 ) \X Ca ) \X Cb ), 

and using manipulations analogous to the lepton-collider case one obtains 


(2.38) 


1 


da 


dr± 2 E'l 


CM 


I j_dT C2 j_d J T Ca j_dT Cf! j_dT s j_S(Tj_ T s _|_ T Ci J_ T c 2 J_ r c a ± T c b ±) 

X ( ^2 ^°l Xc l’ S l^ c l) Xci ,5' |0) 5(t c1 ± - Tj_(X Cl )) 

\ X a 

X | y " (0| Xc 2 ,7 |*c 2 > (X C2 | Xc 2 , 7 ' |0) fi( T C 2 -L ~ r±(X C2 )) 

\W 2 

X ( (-fol Xc a ,a \X Ca ) (X Ca \ Xc a ,a' | Pa) <5(r Ca _l_ ~ T\_(X Ca )) 


J c 


X -w I Xc b ) (X Cb \Xc bl p> |A) <W - T ± (X c J)j C.C} 

X ^35 J/ (r s± ) 7 ^7^7^7^^(2vr) 4 «5( 4 )(P a + P b - p x ) 0(Q± ~ Qo). (2-39) 

In writing the above formula, we have made use of the fact that all four collinear matrix 
elements are color-diagonal and have replaced 


1 


Xc m ,aX Cm ,a' ~^ jy ^ J Xc m ,aXc m ,a' > 


(2.40) 


where it is understood that the color indices of the fields on the right-hand side are con¬ 
tracted. The resulting color contractions then act on the soft Wilson lines and produce the 
soft function 


A7 


Su(r s± ) :=-£ <0| [ylTiY Cb \.. \X„) (X s 

x S(r s± - T ± (X S )), 


YlTjY C2 


J Ik 


YlTjY Ca 


■ Cb 


J0> 


(2.41) 
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which, as anticipated, is a matrix in the color indices I and J . For the (/(/-channel we are 
considering here, there are two possible color structures T\ = t A and T 2 = I, see eq. (2.36), 
and the soft function is therefore a 2 x 2 matrix. 

In a next step, we again expand out the power-suppressed components in the mo¬ 
mentum conservation 5-function. After this step, only the large components of the four 
different collinear fields remain in the 5-function. We already anticipated this step in our 
definition of the soft function in eq. (2.41), which was defined without a constraint on the 
soft momenta other than their thrust. Next, we take the matrix elements of the collinear 
fields along the beam which then define quark beam functions B q / c .±{r Ci ±, Xi), according 
to 


Bq/c a ±(T Ca ±,Xa) := ^ 5(r Ca _L - r±(X Ca ))6 (n a ■ P a (1 - x a ) - n a • p Xca ) 


x c 


W 


x (Pa\ Xc a ,ct(0) I X Ca ) (X Ca \ Xc a ,a' (0) |-Fa) > 

Bq/ Cb ±{r Cb ±,x b ) :=J2^(T Cb ± - T±(X Cb ))6 (n b ■ P b (1 - x b ) - n b ■ p Xc ^ 
x Cb 

x (Pb\ Xcf,,/?(0) I X Cb ) (X Cb \ Xc b ,/3'{0) \P b ) , (2.42) 


where n a ■ P a = n b ■ P b = Ecm- We can introduce the integrations over the momentum 
fractions, x ab , in the expression for the cross section by writing 


1 = ( n a -Pa)(n b -P b ) 


J dx a dx b 5(n a -Pa( 1 


Xa)-n a -Px Ca )S(n b -P b (l-x b )-n b -px Cb )■ (2.43) 


We then have 

~j ~ = 4~/y3 / dr Cl± dr C 2 ± dr Ca± dr Cb± dr s± S(r ± — t s ± — t C i ± — t C2 ± — r Ca _ j_ — t C(j _l) 
x dx a dx b B q/Ca± (r Ca± ,x a )B q/cb± (r Cb± ,x b ) tfg'Ys'y 


X HP <°| Xci,S \X C1 ) {X Cl \x Cl ,5' |0> 5(t c1± - T± (X C1 )) 

\W-1 

X ^ <01 Xc 2 ,7 \Xc 2 ) (X C2 \ Xc 2 , 7 ' 10) 5(t C2 _l - r ± (X C2 )) j CjC*j Sjj(t s± ) 
x (27r) 4 5 (4) (x a P a + x b P b - p Xs - Px ci ~ Px C 2 )0{Q± ~ Qo), (2.44) 


I Tl Tl 1 ' -\-TV nT' 

where gj^ u = - b 2 b a . At this point, the cross section has the form which is usually 

obtained in perturbative QCD, except that the initial state is described by beam functions 
instead of PDFs. To make the following discussion similar to the lepton collider case, it is 
convenient to now perform a boost to the frame where the jets are back-to-back, i.e. the 
frame where 6 \ = 62 . Up to terms suppressed by powers of t±, the total momentum of the 
partons that enter the hard interaction in the original lab frame is given by 

LL -ll 

Ti TV 

q^ = x a Pa + x b P b = X a EcM+ x b E C M-^~, (2.45) 
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where we used that rib = n a . At leading power a boost along the z-direction is thus 
sufficient to make the jets back-to-back. We need to perform a boost such that the total 
momentum of the partons that enter the hard interaction is given by 

Li -11 

yi' 71 

t ^ r = Q f + Q y + 0(Qt ± ), (2.46) 

where we denote quantities in the boosted frame with a hat. The plus and minus compo¬ 
nents of the c a particles transform under a boost in the z direction according to n a ■ p —> 
e^n a ■ p, n a ■ p —>• e~^n a ■ p, therefore the boost we need is given by 

£ = - In — ; Q = EcMy/XaXb- (2.47) 

2 x a 

Since the constraint for transverse thrust only involves momenta transverse to the beam 
direction, the expression for the cross section in eq. (2.25) is invariant under the boost. 
The beam functions are invariant and can be evaluated in the lab frame, while we will 
adopt the boosted frame for the jet and soft functions. From now on we will exclusively 
work in the boosted frame, and for simplicity drop the hats on the momenta. We then 
rewrite the momentum conservation ^-function as in eq. (2.16) and, also as in the lepton 
case, introduce an integration over the angle 8 of the jet in the boosted frame by writing 

1 = Jd 3 n6^ 3 \n — ni) = (2ir) j dcosd ^ 2 \px C i ) ■ (2.48) 


Finally, we take the collinear matrix elements, which are the jet functions defined in 
eqs. (2.19)-(2.20), where Q is now given by eq. (2.47). Putting all the ingredients together 
we obtain 

7 /» 

• / dl~ Cl j_ dT C 2 J_ dr s± d,T Ca j_ dT Cb j_$( tj_ T s j_ 7"c2-L ^~C1_L ^~Cb- L ^”c a -L) 

X I dXa I dXb ^ 9 / Ca - L ^ Tca - L ’ Xa )Bq/ c b±(. T Cb-Li Xb ) Jc 2l( Tc 2-L) Jc i-L( T ci±) 


dr I d cos 8 


x Hij(Q, 0) Sji±(t s _l) 9(Q± - Qo) 


(2.49) 


where we defined the hard function as 

Hu(Q,0) = ^^-(i + cos 2 8)C I C*J. (2.50) 

Note that, unlike in the lepton-collider case above, we have not normalized it such that the 
entries of the matrix are only ones and zeros at tree level. The angle 8, which corresponds 
to the angle between jet 1 and the beam in the partonic center-of-mass frame, is typically 
not measured experimentally, but it is convenient to keep the expressions differential in 
cos 8 for the present discussion. 

We can verify that we reproduce the lowest-order cross section for two-jet production 
by evaluating the ingredients at leading order. The LO hard coefficient is Cj = g 2 s (1,0 )/Q 2 ■ 
Its denominator Q 2 is from the hard gluon propagator which is integrated out when match¬ 
ing onto SCET. The relevant entry of the soft-function matrix is 5m(r s i) = fd(T s ±). 
The LO jet functions are equal to <5(r Ci j_) and the beam functions are given by the standard 
PDFs multiplied by the same (5-function. Putting the ingredients together, and integrating 
over Tj_, one reproduces the standard LO partonic cross section for qq —> q'q'. 
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2.3 Collinear anomaly 

For leptonic collisions, the theorem eq. (2.22) achieves the complete factorization of scales 
we are aiming for, but the same is not true for the hadron collider formula eq. (2.49). 
The reason is that the beam and soft functions, as given above, are not individually well 
defined. Beyond leading order, one encounters phase-space integrals that are ill-defined 
within dimensional regularization. To properly define them, one needs to introduce an 
additional regulator. The convolution of the regularized beam and soft functions is then 
regulator independent, as it needs to be since the cross section is a physical observable, 
but contains an additional dependence on the hard scale, on top of the one encoded in 
the hard function. This effect, i.e. the appearance of a hidden additional dependence on 
the hard scale, is called the collinear anomaly [13]. One encounters it when the effective 
theory contains collinear and soft modes which have virtualities that are parametrically 
of the same order. Consistency conditions restrict the form of this additional hard-scale 
dependence to be a pure power to all orders in perturbation theory [13, 21]. The additional 
dependence can be factorized and in this way one achieves the desired separation of scales. 
One can show that for transverse-momentum dependent quantities, such as the one we 
consider here, only real-emission diagrams need additional regularization and one can use 
a simple analytic regulator to render them well defined [22], In this way, gauge invariance 
and the eikonal structure of the soft and collinear emissions in the effective theory are 
explicitly maintained. Alternatively, one can regularize the soft and collinear Wilson lines 

[23]- 

To discuss the form of the collinear anomaly it is convenient to perform a Laplace 
transform of the cross section, since the Laplace-transformed expression is a simple product, 
rather than a convolution, of the jet, beam, and soft functions. We write 

da 

dr±d(cos 9)dx a dxb 

= Hu(Q , 9) S ji± (k, 9) j' C2 _l(/v) j ci ± (k) B q/Ca _ l (k, x a ) B q / Cb± (K, x b ), (2.51) 

where z = l/(e lE n). The factor e 1E is included to avoid a proliferation of such factors in 
the Laplace transforms of the jet, soft, and beam functions, which are indicated by a tilde. 
For example, for the jet function we define 

jc 1 ±(«) = J dr±e~ T±z J Ci ±(t±), (2.52) 

and analogously for the other ingredients. For later convenience we also define the dimen¬ 
sionful variable R = 2k Q sin0, which is of the order of the soft energy scale. As mentioned 
before, the soft and beam functions in eq. (2.51) depend on the analytic regulator, and 
contain hidden dependence on the hard scale. Their product is regulator independent and, 
to all orders in perturbation theory, has the form 

~ ~ ~ / q2 

Sji±(k, 9) B q/Ca± ( K , x a ) B q/Cb± (n , x b ) = Wji(k, 9, x a , x b ) , (2.53) 

\ c o K / 
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where 


Wj^K^e^Xa^Xb) = Sj I± (k, 6) Bq/^i^Xa) Bq/ Cb± (K,X b ). (2.54) 

All the dependence on the hard scale Q 2 is now manifest. The structure we find here closely 
resembles the one encountered in previous applications involving the collinear anomaly, see 
e.g. refs. [24, 25], and the derivation of eq. (2.53) parallels the steps presented there. 
We call the function F±(k) the anomaly exponent, and Wjj(k) the remainder function. 
We split the remainder function into three parts in eq. (2.54). This decomposition is not 
unique, but it is useful because it shows that the dependencies on x a , x b and 6 factorize. 
The constant Co is conventional and will be given in section 3. We can now write the final 
factorized form for the Laplace transformed cross section. For a general partonic channel, 
ab —> 12, the corresponding factorization formula reads 


t ab ^ 12 (n) = 6) 


ab—>12 1 


Q' 2 

cZr 2 


2 x -Ff («0 


s" 2 (M) 

B a /c a ±_ («, Xa) B b / Cb _ L (k, x b )j Cl _L (k) j C2 ± (k) , (2.55) 


where we explicitly indicated the partons upon which each term depends. Equation (2.55) 
is the main result of this work. In section 3, we compute all the ingredients that enter in 
this formula at the one-loop level. 

As discussed above, the cross section must be independent of the analytic regulator, 
while the individual pieces have divergences as the regulator goes to zero. This cancellation 
of divergences is non-trivial in our case because the soft function is a color matrix, while the 
beam functions are color-diagonal. It is interesting to look at the structure of the analytic 
divergences in the soft function in detail to check how they cancel the divergences in the 
beam functions. In our computations in the next section, we will use the standard form of 
the phase-space regulator [22] 


Jd d kS(k 2 ) 9(k°) 


J d d kd{k 2 )0(k°) 


V 

rib * fc 


a 


(2.56) 


with n b ■ k = ko + k z , but for the discussion in this paragraph, it is convenient to introduce 
the analytic regulator as follows [26] 


J d 4 kd(k 2 ) 6(k°) 


V 

rib * & 


a 

0(n b ■ k - n a ■ k) + 


v 

n a ■ k 


0 

0{n a ■ k — n b ■ k) 


(2.57) 


where h a = n b and n b = n a - This second form distinguishes the divergences from left- and 
right-moving particles and is symmetric if one chooses a = j3. The implications of having 
separate regulators for left and right sectors were first discussed in ref. [27]. Since the 
components k+ and scale differently in the collinear regions, the regulator (2.57) must 
be expanded in the small components. In the collinear region c a we have n b ■ k n a ■ k 
and the term involving the regulator /3 vanishes upon expanding the ^-functions, while the 
a term vanishes in the collinear region c b . We conclude that the divergences in the c a 
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beam function only involve the regulator a and vice versa for the c b beam function. The 
question is then how the cancellation can be possible, given that the soft function depends 
on the color charges of all four particles participating in the hard scattering, while the 
beam functions only involve the color generators of the incoming partons. 

To see how the cancellation arises, it is instructive to consider the one-loop case, where 
the soft function is given by a sum of diagrams in which a gluon is emitted from leg i and 
absorbed at leg j (see figure 7 below). The amplitude squared associated with this process 
has the form 

< 2 - 58 ’ 

We use the notation of the color-space formalism [28], in which T) denotes the color gener¬ 
ator associated with parton i. Only the contributions which involve leg a (i = a or j = a), 
can give rise to divergences in the regulator a. Computing the corresponding diagrams 
explicitly, we find that they have the structure 

- (T a • T b + T a • Ti + T a ■ T 2 ) = - (-T a 2 ) , (2.59) 

a a 

where we have used color-conservation JT T) = 0. The quantity T 2 = C a is the quadratic 
Casimir operator of the representation of parton a. Similarly, we find 

\ (T a ■T b + T b .T 2 + T b -T l ) = -p (-T 2 ) . (2.60) 

From eqs. (2.59)-(2.60) one clearly sees that the divergences in each hemisphere can cancel 
with the corresponding analytic-regulator divergences in the beam functions, which are 
proportional to the corresponding Casimir. The same structure must also arise at higher 
orders. Factorization thus imposes nontrivial constraints on the structure of the divergences 
in the analytic regulator, similar to the constraints it imposes on the infrared structure 
of the scattering amplitudes [20, 29-31]. The above structure implies that the collinear 
anomaly has the form 

F?(k)= {Ca+ 2 Cb) F ± (K). (2.61) 

The anomaly is thus given by a universal function Fj_(k). Casimir scaling of the anomaly 
was observed earlier in the simpler case of transverse-momentum resummation [13] . 

3 One-loop ingredients 

In this section we compute all the different ingredients that appear in the factorization for¬ 
mula at the one-loop level. We work in dimensional regularization with d = 4 — 2e dimen¬ 
sions, and a s = (7 2 /(47t) is always understood to be the coupling constant in the MS scheme 
at scale fi. The relation to the bare coupling is = Z a a s fj ? £ , with fi 2 = / u 2 e 7E (47r) _1 , 
and Z a = 1 at the order we are working. We expand the anomalous dimensions appearing 
in our expressions according to 

tx = E7»,(£)” + ‘ ; 7cu,p = f;r„(^)" +1 . (3.1) 

71—0 71—0 

where the coefficients T n of the cusp anomalous dimension are collected in appendix A. 
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Figure 5. Next-to-leading order real-emission diagrams for the quark jet function (first line), and 
for the gluon jet function (second and third lines). The red vertical lines represent the final-state 
cut. The crosses indicate the collinear Wilson lines. 


3.1 Jet functions 

For simplicity we denote n Ci —> n and Cj — > c in this section, since the diagrams only involve 
one collinear direction, and represent the quark jet function by J q ±(T c ±) and the gluon jet 
function by J 3 _l(t c _l). The definition of the quark jet functions was given in eqs. (2.19) and 
(2.20). For hadronic collisions, we also need the gluon jet function whose definition reads 

xab 

2Q{27T)d -i (-&) Jg±(r c± ) := £ <0| A a c t m (0) \X c ) (X c \ A b cTu ( 0) |0) 

x c 

x 6(t c _ l (X c ) - t c ±)5(Q - ft ■ Px c )S {d ~ 2) (p Xc ,T), (3.2) 

where Aj^ is the SCET gluon jet field (normalized such that Aj/i = A-j^ at lowest order 
in the coupling constant), and gj w = — n ^ n "'+ n ^ n ^ ^j xhe normalization is chosen such 

that Jg±(r c ±) = 5(t c ±) at lowest order. To ensure that the same is true for the quark jet 
function in d dimensions one has to replace (27r) 3 — > in the definitions in eqs. (2.19) 

and (2.20) and work with a (d — 2)-dimensional transverse-momentum h-function. 

The diagrams that contribute to the jet functions at the one-loop level are shown in 
figure 5. We only show real-emission diagrams since the one-loop virtual diagrams are 
scaleless and vanish in dimensional regularization. It is convenient to follow ref. [35] and 
perform the calculation of the jet functions in light-cone gauge n ■ A(x) = 0. In this gauge, 
the collinear Wilson lines which multiply the fields Xc an d Ac become trivial, therefore 
only the first diagram contributes to the quark jet function and only the first two diagrams 
in the second line to the gluon jet function. The relevant phase-space integrals can easily 
be performed in d-dimensions and we then obtain the bare quark jet function at one loop 
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as 


J 


bare 

<?-L 


^ w_ ^ ^ (4- £ )2- 3 + 2 WT(2- £ ) 1 fr c± Q^si^6\ £ o , 

( c±) - ( c±) s F e r(i- e )r(|- e ) t c _l \ ^ ) ' (3 ' 3) 


The explicit expressions for the two diagrams contributing to the gluon jet function are 
given in appendix B. 

The bare jet function has divergences for e —> 0. Since the jet functions are distribu¬ 
tions in t c _|_, it is convenient to perform a Laplace transform 


ji±(L,n)= / dr c _ L e Tc±/(KeJE ' 1 Ji±(r c _ j_,/t), 

Jo 


(3.4) 


where, for later convenience, we have written the Laplace transform as a function of the log¬ 
arithm L = ln( 4fc< ^ 2 * n 9 ) °f the Laplace variable k. The renormalized Laplace-transformed 
jet function is related to the bare one according to ji± = Zj t ± j l b j) rc . Since the bare func¬ 
tion is //-independent, the renormalized jet function, jjj_, and Zj t ±_ fulfill both the same 
RG equation 

(3.5) 


d 


— 2 Ci 7 cusp L ji±(Li fJ-) 


din fi“ 

where 7 CUS p is the cusp anomalous dimension, which gets multiplied by the relevant Casimir 
operator C q = Cp and C g = Ca- Solving the RG equation for Zj t ± one gets 


1 7 ° is 

in Zf i = — 
4 t r 


+ -(Ci r o L + 7q1) 


(3.6) 


Expanding and Z 3i ±_ in e and a s , we find that the divergences in ji± cancel with 


7oi=7o 9 = -3 C F , 


7 0 1 = 7o = -A) • 


(3.7) 


The one-loop anomalous dimensions are thus identical to the inclusive jet-function with 
our choice of the scale of the logarithm L. The renormalized function ji± is given by 


ji±(L, At) — 1 + ^ (CiT q— + 7 q‘ ± L + c(' ± 


with 


a± = Cf[7-^- 


3.2 Beam functions 


67 47t 2 


'■J- = Ca I T “ T j “ ¥ Tf "/' 


20 , 


(3.8) 


(3.9) 


The definition of the quark beam function was given in eq. (2.42) in section 2.2. In addition, 
we also need the gluon beam function which is defined as 

Bg /caX _(TcaX,x a ) := -x a (n a ■ P a ) £ (P a \ A^ T { 0 ) \X Ca ) (X Ca \ A a CaT ^0) \P a ) 

x 5(r Ca ±(X Ca ) - T Ca ±)S ((h a • Pa)(l - Xa) - n a ■ Px Ca ) • (3.10) 
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Figure 6. One-loop diagrams contributing to the matching coefficients Z q <- q (top row), I g< _ g 
(middle row), and Z q <- g (bottom row left) and X g ^ q (bottom row right) of the beam functions. The 
red vertical lines represent the final-state cut. 


For perturbative values of tj_, we can compute the final-state radiation in perturbation 
theory and match the beam functions onto standard PDFs, fj/N ( x ) for the hadron N, 
convolved with a perturbative coefficient, t±). The corresponding matching relation 

takes the form [13, 17, 36] 

/•l dx 

&i/c k . l( t ±! 0 = £ / (£/x,t±) fj/N k (x). (3.11) 

3 ^ X 

At lowest order in perturbation theory the matching coefficients take the form 

%i<-j (x , Tj_) = 6ijS( 1 - x)6(t_ l ) , (3.12) 

so that the beam function reduces to the PDF times 5(r j_). When computing the matching 
coefficients, it is convenient to evaluate the beam-function matrix elements eqs. (2.42) and 
(3.10) with quark and gluon states instead of hadrons. In this case the PDFs are trivial 
fj/k( x ) = djkS(l — x) and the partonic computation directly yields the matching coefficients 
2^ fe (£/x,T]_). 

As explained in section 2.2, the transverse-thrust beam functions are not well defined 
within dimensional regularization and require an additional regulator. We regularize them 
using the analytic regulator of ref. [22] which modifies the phase-space integration measure 
as follows 

J d d k5(k 2 ) 9(k°) -> Jd d k 5(k 2 ) 9(k°) ^ . (3.13) 

Note that the beam functions in both c a and q, directions are regularized with the same 
regulator. To discuss the renormalization of the beam functions, and the cancellation of 
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the dependencies on the analytic regulator, we again perform a Laplace transform 


Bq/c a ±(K,X a ) = / d,Tj_e T ^ /{K ^ E) B q/Ca± {T L ,Xa). 


(3.14) 


We have discussed the cancellation of the divergences in the analytic regulator in sec¬ 
tion 2.3. The divergences of the soft function with the form of the analytic regulator given 
in eq. (3.13) can be obtained from those in eqs. (2.59)-(2.60) by replacing (3 —a. We see 

that the analytic divergences cancel in the soft function for channels involving the same 
color representation of the incoming partons, such as qq and gg. In these cases, the remain¬ 
ing divergences must cancel out in the product of the Laplace transformed beam functions 
which takes the form 


Bq/ Ca ±(K,x a ) B q/cb± {K,x b ) = 


Q~ 

cln 2 


2 \ ~CfF±{k,h) 


B q/c a ±(^X a ,g) Bq/ Cb± (K,X b ,n) , (3.15) 


and analogously for the gg channel. Computing the product of the Laplace transforms of 
the quark beam functions, we find that this is indeed the case, and obtain the corresponding 
anomaly exponent and refactorized matching coefficients. The diagrams that contribute to 
it at one loop are shown in figure 6. The diagrams need to be evaluated both in the c a and 
Cb sectors because the analytic regularization is not symmetric, as we stressed above. The 
explicit expressions for each of the individual diagrams are listed in appendix B. Proceeding 
as discussed, we get 

rv _ 

(3.16) 


Ole 


f±(k,h) = — r 0 T _ L , 


with L i = In 4^2, with lnco = and where G ~ 0.915966 is the Catalan constant. The 
characteristic scale of the logarithm is R = 2 kQ sin#. The refactorized matching coefficients 
are defined as 

rl d z 

Bi/ Ck ±(r±,£,lJ) = V / — Ii^j(€/z,T ± ,iJ,) fj(z,n) , (3.17) 

Je z 

in analogy with the corresponding matching coefficients and are obtained by evaluating the 
beam functions with partonic states. At one-loop order, they have the general form 


Ii<-j(x, k, /j,) = 8(l-x)6 , 


y 


a.s 


Li 


1 + t 


a s 

+ ~j— 

47T 


-v^lj(x) I ^+n i ^ j {x) 


and contain the Altarelli-Parisi kernels 

1+ x 2 


vW q (x)=4 C F {^~ 


— x 


V$L g {x) = 4 T f (x 2 + (1 - x) 2 ) , 

?>«(,) = 8 C A [ X 1 ~ X 


+ 


.(! - x)+ 

pW (x)=4C F 1 + (1 y l)2 , 


+ x(l — x) 


+ 2/3 0 d(l - x ), 


(3.18) 


(3.19) 
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which involve plus-distributions, and remainder functions 


n q ^ q (x) = -C F + 25(1 -x) + 2C f (1 - x ), 

7Z q ^g(x) = 47>x(l - x ), 

Kg^g(x) = -C A ( y + + 2i? ) -K 1 - *), 

lZg^ q {x) = 2 C F x . (3.20) 

The numerical value of the constant F in the equation above is F ~ 8.20629, an analytic 
expression in terms of Lerch’s <f> function is given in the appendix. 

The anomalous dimensions of the anomaly exponent and the refactorized beam func¬ 
tions are given by 


din fi 


F± = 27 , 


cusp 5 


din /a 


-B 


i/ck -L 


r 2 K 2 

l < '0™ o rs 

-CB Z 7cusp m —^2 - 2 Ti 


Bi 


B 


i/c k ± 


with 


c Bi = Ci and 7 ^ = 7 ^, 
according to our one-loop results above. 


(3.21) 


(3.22) 


3.3 Soft functions 

The soft function is, in general, a matrix in color space defined in the color basis adopted 
for the hard function. To discuss the color structure independently of the basis choice, it 
is most convenient to use the color-space formalism [28, 37], where the hard function is 
written as 

\C) = J2Ci\ I ), (3.23) 

1 

with |J) a basis of color states. For the qq q'q' channel discussed above, the basis 
contains the color-singlet and the color-octet operators defined in eq. (2.36). In the chosen 
basis, the soft function acts as follows 

{C\S(T a± )\C) = C} (I\S(t s± )\J) Cj = Cj S IJ (r s ±) Cj = H n S u ( t s± ) . (3.24) 

We can write it as 

Su(r s± ) = ±r J2 <0| Wj \X 8 ) (X 8 \Wj |0> 5(r s± - t ± (X 8 )), (3.25) 

c x B 

where Wj are combinations of Wilson lines, in the adequate representation for the partonic 
channel we are studying. The explicit expression relevant for the qq —> q'q' channel was 
given in Eq. (2.41). 
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Figure 7. Next-to-leading order real-emission diagrams for the soft function. The red vertical lines 
represent the final-state cut. Mirror diagrams and those that vanish because the gluon attach to 
two Wilson lines in the same direction are not shown. In the lepton-collider case the a and b Wilson 
lines are absent, and only the last diagram is possible. 


The diagrams that contribute at the one-loop level are shown in figure 7. Once again 
only real-emission diagrams can give non-vanishing contributions in dimensional regular¬ 
ization. The soft function up to one-loop order is therefore given as 

Sjj(t) = 8{t)D)Y + Y, (3.26) 

ij 

where i,j = a, b , 1, 2, with a and b the partons in the initial state, and 1 and 2 the final-state 
jets. The color structures at tree level and one-loop order are 

Dt/r = W C {I \ J ) ’ and D IJ = ^(I\(-T r T J )\J). (3.27) 

The Djj encode the color factors, are matrices in color space, and depend on the partonic 
channel we are considering. For the qq -type channels they are 2x2 matrices, for the qg -type 
channels 3x3 matrices, and for the gg channels 9x9 matrices. Explicit expressions for all of 
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them were given in ref. [38]. 3 For leptonic collisions the color factors are numbers because 
the process only involves two color-charged legs which fulfill — T\ ■ = T/ = Tf = Cf- 

The soft integrals I %3 only depend on the directions to which the soft gluon in figure 7 
attaches, and not on the representation of each Wilson line. They are given by 


lij — 


drk 


I T^-^s(k 2 )9(k u ) , n ;' nj ,, 

J (27T) d — 1 V ; V J (m ■ k){nj ■ k) 


8 t s _l - 


1 


Q| sin 9 1 


— |nj_ • k _l 


(3.28) 


with 


K = (!, o, 0,1), = (1,0, 0,-1), < = (l,n), n£ = (l,-n), n = (sin0, 0, cos0). 


For leptonic collisions, only I \2 arises, since there are no colored particles in the initial 
state. The Uj integrals are in general not well defined within dimensional regularization, 
and an additional regulator is required. In order for the divergences to cancel, we need to 
use the same analytic phase-space regulator as for the beam functions. The regularized 
versions of the lij integrals have an additional (v/{kQ + k z )) a factor in the integrand. 

To compute the integrals lij, we find it useful to perform a change of variables 
( k ±, k z ) ->• ( x,y) with 


k _l = 


t s± Q\ sin6>| 


_ r s _L(5|sin0| n 

"'Z — V- 


(3.29) 


The variable x G [0,1], while y 6 (—oo, +oo). To extract the divergences in a , it is 
convenient to introduce the variable u = 1/(1 — y + y/\ + y 2 ) G [0,1]. Parameterizing the 
integrals in this way, it is more or less straightforward to extract the divergences and to 
compute the remaining finite integrals numerically, but we have not managed to evaluate 
the finite parts analytically. The nontrivial soft integrals can be written in the form 

I vi = 


2/„i = 

2/&2 = 


2 

n a ■ ni 


2 N(l + +1_+I' + -1 , _), 

7T~N (1+ - 1_ + l’+ + 1'J) . 


(3.30) 


Below, we will discuss the evaluation of the integral Zq in detail. The normalization factor 
N and explicit expressions for the integrals Z±, Z± are given in appendix B. The three 
remaining soft integrals are given by 


lab 0 , Ia2 — y-rc—Q i ^bl — nr—6 

The integral I a b is zero because it is proportional to the scaleless integration 

(y+ + y 2 ) 


f 


hi 


dy- 


\/l + 


y z 


= / dutr 1 -a (i-u) 
Jo 


— 1 +Q _ 


= 0. 


(3.31) 


(3.32) 


li Note that the matrices in ref. [38] are defined without the l/iV c prefactor in eq. (3.27). Some of the 
entries in the matrices for the gg channel in this reference are interchanged: the entries 4 and 19 (20 and 
65) need to be swapped in the tree-level (one-loop) matrix. We thank Jan Piclum for pointing this out. 
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The I \2 integral is well defined without the additional analytic regulator, since it does 
not involve the directions a and b of the initial-state partons, which involve collinear modes 
of the same virtuality as the soft modes. The soft function for leptonic collisions is there¬ 
fore well defined within dimensional regularization, and there is no collinear factorization 
anomaly, as expected from the arguments in section 2.3. To illustrate how the soft integrals 
are computed, we explain now how to obtain the result for Xq, which is given by 


Xq = / dx (2 — x) 


2 £ x 2 + £ 


-ff y/l + y 2 1 + y 2 - (ycos<9 + (1 - x) sin#) 2 ' 


To arrive at the above form from eq. (3.28) one writes the integration measure as 



dko dk z dk± k d ± 3 d(cos ip) (sin Lp) d 5 fld -3 


d —5 i 


(3.34) 


and uses 5(k 2 ) to perform the k° integration and the transverse-thrust constraint to obtain 
the integral over the angle <p. Finally, one performs the change of variables (k ±, k z ) —> ( x , y) 
specified above. The next step is the y integral, which does not depend on e. To compute 
it, we write 

poo 

/ dy- 


where 


1 


1 


/ 


— OO 

oo 


y/l + y 2 1 + y 2 — {y cos 9 + (1 — x) sin 0)* 


dy 


1 


1 


i/l + y 2 (yi ~ V 2 ) sin 2 9 \y -yi y-y 2 


l 


l 


dy 


2 i 


Im 


V 7 ! + y 2 ( 2/1 - V 2 ) sin 2 9 \y-yi 


(3.35) 


„ (2-x)2x2 ,(2-x)2x2 

y i = (1 — x) cot 9 + i- - ; — - ; yi = (1 — x) cot 9 — i- 


sin0 


sin0 


(3.36) 


are the roots of 1 + y 2 — (ycos9 + (1 — x) sin0) 2 ; note that y\ — yi is purely imaginary for 
0 < x < 1. Using 


dy 


1 


1 


1 


log 


'A-+A 2 


we find 


Xo = (sin 9) / dx (2 — x) £ x 

Jo 


sJY+tfy- A v / TT^42 6 yA + y/l + tf) ’ 

yi - V 1 + y'i 


(3.37) 


— 1 — —1+£ 


Im 


1 


: In 


\/l + y i \ 2/1 + \A + '” 2 


y 


(3.38) 


The divergences in the integral over x only come from the x 1+£ term, and the result for 
Xq as a series in e can be easily found by using the expansion 


l+£ 


= itf(s) + 


+ £ 


lnx 


x 


+ 0(s 2 ), 


(3.39) 


J + 
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where the plus distributions are defined by 

+ g(x)dx= [ f(x){g{x) - g(0))dx. (3.40) 

Jo 

The integrals 1+ and X' + are finite as a —> 0 and can be evaluated in exactly the same way 
as Iq. A divergence in the analytic regulator a arises only in the integrals and I'_ in 
eq. (3.30), which are simpler than the plus-type integrals so that the y-integration can be 
performed analytically for arbitrary a. The divergent terms of the integrals in eq. (3.30) 
are given by 



h 2 = 

Ota 

1 

7 r 

T s ± 

24 1 = 

Ot 8 

1 

7 r 

r s ± 


X 1 

W 

27f,2 = 

Ota 

1 

7T 

r s ± 


x! 

(1 _ 


/ 4Qt s _l 


* o \ —2s / 1 
sin 6 \ / 1 


— 2 In sin # + 0{e) ] , 


/ 4Qt s _l 


sin# 


—2e 


9 ( 2Qt s ± sin 9 

2 + V v 


2 ( — + H ) + 0(e, a) ] , 


/ 4Qt s± 


sin# 


—2e 


- - ( Q Ts± 


sin# 


2 V 2z/ 


2 ( —h H ) + 0(e , a) 


(3.41) 


where the numerical value of the constant H is H « —1.85939. We refer to appendix B for 
the full expressions, including the finite parts. 

Due to the structure of the Du color matrices, see e.g. ref. [38], the I- L j integrals always 
appear in the combinations I \2 + I a b , lai + h 2 , and I a 2 + hi for the qq- and ^-channels. 
As we can see from eq. (3.30), and explicitly in eq. (3.41), in these combinations the 1/a 
divergent terms cancel, as expected from the discussions in section 2.3, and the explicit 
calculation of the beam functions in the previous section. For the qg- type channels, on the 
other hand, the structure of the Du matrices is such that some individual hj integrals 
appear, and not only the above combinations. Therefore, in the qg- type channels, the 
divergences in a do not cancel within the soft function itself, but one needs to combine 
it with the beam functions to obtain a regulator independent result, again as discussed in 
section 2.3. 

The RG equation of the soft function is given by 

s a ^— s a ^( r a6 ——>i2t\ c a b —^12 /q 421 

d ln/i IJ± ~ IL± V s )lj~\ s Jim mj± ’ 1 ' ’ 

with 

(r^ 12 ) LJ = f 7 cusp cs In AkQs ™ 2 9 + Su + 7cusp Mu(0) ■ (3.43) 

The coefficients eg, 7s j_, and the remainder function Mjj(9) depend on the partonic chan¬ 
nel. Up to two-loop order the remainder function is fixed by general constraints on the 
structure of soft anomalous dimensions, see section 4. From our results above we get 


eg = —(Ci -I- C 2 ) ; 7 os_L = 0. 


(3.44) 
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3.4 Hard functions 


Like the soft function, the hard function Hjj in eq. (2.50) is a matrix in color space. Up to 


the conventional normalization factor, it is given by the Fourier transform of the matching 
coefficients Cj of the QCD currents to the SCET operators, according to 


Hu oc CjCj . 


(3.45) 


The hard function is extracted from the results for the amplitudes for 2 —> 2 processes in 
QCD. This was done at the one-loop level in ref. [18], where the RG evolution of the hard 
function was also considered. Recently, the two-loop result for the hard function has been 
presented in ref. [19], based on the results for the QCD amplitudes of refs. [39-43]. The 
results of ref. [19] are conveniently given in an electronic form, and can be readily used. 
We refer to refs. [18, 19] for the explicit results for the hard functions, and do not copy 
them here. 

4 Resummation 

4.1 Renormalization group equations and scale independence 

The cross section, or equivalently its Laplace transform, must be RG invariant, i.e. it must 
be independent of /x. That this is indeed the case at order a s can be explicitly verified 
with the one-loop results presented above. The Laplace transform of the cross section was 
given in eq. (2.55). For concreteness, we now consider the cross section in the qq channel 
whose Laplace transform is given by 



(4.1) 


where Hu := H qq ~* qq and Sij := S qq ~^ qq . The cross section must satisfy 


dt 


(4.2) 


din fi 


which implies relations among the anomalous dimensions of the different ingredients. Their 
RG equations can be written as 



(4.3) 
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The general structure of the hard-function anomalous dimension was derived in [20, 29- 
31]. Following [18], we rewrite the anomalous dimensions as a diagonal contribution and a 
remainder Mi j (cos 6) 

(r H )u = (jcusp^ In $IJ + 7cus P Af/j(cos6'). (4.4) 

The remainder depends on ratios of Mandelstam variables, which can be rewritten in terms 
of the scattering angle 6. In order for the related angular dependence to cancel, the soft 
anomalous dimension must involve the same remainder Mij (cos 6) 


(r s)u 


^TcuspdS In 


4 kQ sin 2 6 


+ 7s± 


Sij + 7cus P A4/j(cos 9). 


(4.5) 


For the case of the soft function, the angle dependence arises via scalar products of the 
light-like reference vectors defining the Wilson lines. In addition, eq. (4.2) also imposes 
some constraints on the coefficients and anomalous dimensions appearing in the diagonal 
parts of the equations above, which are 


c H -cs-2c Bq -2c jq =0 ; lH - ls± - 2 7 f* - 2 7 J = 0. (4.6) 


These conditions for the diagonal and non-diagonal parts are verified by our one-loop results 
in the previous section, which provides a check of the computations. To verify the scale 
independence, one uses relations such as 

u _ rig ■ ri 2 _ s (~t) _ n\ • n 2 n a • m 1 - cos 0 

t n a -ni C ° 2 ’ u 2 (n a -n 2 ) 2 (l + cos 0 ) 2 ’ 

where s = (p a +Pb) 2 , t = {p a — pi) 2 , and u = (pb — Pi] 2 are the Mandelstam variables. 


4.2 Resummation of large logarithms 


With the RG equations at hand, we can now derive general resummed expressions for the 
cross section. To do so, we solve the equations in Laplace space and then invert the results 
back to momentum space, using the technique of [44], All RG-equations in eq. (4.3) are of 
the form 


d 

din p 






Cf 7cusp In 


(V 


+ 7/ 






(4.8) 


where Af is the characteristic scale of the given function, Cf the relevant combination 
of Casimir operators, and 7 / its anomalous dimension. For the hard function one has 
A h = Q, for the jet function A 2 = 4KQ 2 sin 2 0, the soft scale is A s = 4K,Qsm 2 d, and the 
beam functions depend on A& = cqk = 2coft(Jsin0. The RG equations for the hadron- 
collider soft and hard functions are matrix valued; to bring them to the form shown in 
eq. (4.8), one first has to diagonalize the anomalous dimension. 

The solution of the template RG equation (4.8) reads 




exp[— CfS(p f ,p) 


-Vw.ri] (^) • 


(4.9) 
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The evolution factors and Af(nf,[j,) are given in appendix A and the exponent 

rjf = CfA lcusp (jif, /i). To obtain the solution of the RG equation in momentum space, 
one makes use of the fact that at any order in perturbation theory is just a 

polynomial in the logarithm so that one can rewrite 



(V 

di 



di 


Vf 


f(o,,„»,) yyy . 


(4.10) 


After this, the A f dependence, and therefore also the n dependence, is a pure power for each 
ingredient of the factorization formula. Since the cross section factorizes into a product in 
Laplace space, the entire cross section is proportional to a power of k and using the fact 


that 


dr _l e 


T a - l P~ a lE 
-Tj_/(K.e<E) \l_ e 


(4.11) 


j o r (a) 5 

we can then invert the Laplace transform and obtain the resummed result in momentum 
space. Using the above template, we can solve the RG equation for all the ingredients 
in the factorization theorem. By evaluating each one at its characteristic scale and then 
combining them at a common scale //, one resums the large logarithms. Below, we choose 
/i = i-i s for simplicity in order to avoid evolving the soft function. 

The factorization formula in eq. (2.22) for leptonic collisions resembles the cross section 
for ordinary thrust. Its Laplace transform reads 


i(/c) = / dr±e 

Jo 


—rx/(ree 7 -E) 


da 


dr±d(cos 9) 


irN c Q 2 a 2 

2 Q 2 


(1 + cos 2 9)H(Q 2 ,h) 


. 4kQ 2 sin 2 6 \~ f 4 kQ 2 sin 2 6 \ _ / 4 kQ sin 2 6 

x J c ± ( In- 15 -, H Jcl. In-^- ,n s± In-, n 






d 


(4.12) 


where we denote both jet functions by j c ±. We wrote the Laplace-transformed jet and 
soft functions as functions of the logarithm of the arguments so that we can directly use 
the template eq. (4.8) to solve the corresponding RG equations and to invert the Laplace 
transform. After a few simplifications, one then obtains the resummed result 


1 da 
(Tq dr±d(cos 9) 


= g (1 + cos # 2 ) exp [4 C F S(n h ,Hj) - 2A lH (fx h , fx s )\ H(Q 2 ,fi 2 h ) 


x exp 


4 CfS(h s , fij) + 4 A^, (dji ds ) 


3c± y n> -yr + d v-d 3 


91 

dl 


2 \ ‘ZCpA 7cusp (^h J f^j ) 


S' ± (5 f? ,/r s )(4sin 2 0)"- ^ 

T± V ds ) r(r/) 

(4.13) 


where /x/j, /ij and n s are the hard, jet, and soft matching scales, at which the respective 
functions are evaluated. Up to the additional angle-dependence, the result has the same 
form as the one for ordinary thrust derived in [32]. The Born-level cross section is (To = 
—Aj/— and Tj = 4C'p’A 7cusp (/Xj, /j, s ). The factor 3/8 on the r.h.s. accounts for the integral 
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over cos 6 


/ dcosO (l + cos# 2 ) = 

I -1 3 


(4.14) 


At NLL accuracy, one can perform the full angular integral in eq. (4.13) analytically using 
the identity 

r 1 „ . Ptt (1 4- nl 

(4.15) 


? m, _ \/7r T(1 + rj) 


/ d cos 9 (sin^ 6) 71 = - „ . 

1-1 rd + r?) 


To obtain the resummed result also for the hadron-collider case, we start from eq. (2.55) 
and evolve the hard and jet functions from the scales fXh and [Xj down to the scale [x at 
which the soft and beam functions are evaluated. The resummed expression for the hard 
functions in the different channels were given in ref. [18] and we will not reproduce them 
here. The solution for the jet functions is the same as in the lepton collider case. What 
remains is the anomaly times the remainder function Wj b f ^ 12 , given by the product of 
beam and soft functions, see eq. (2.54). To be able to use eq. (4.11) to invert the Laplace 
transform, we follow ref. [34] and write 


Wj b i^ 12 (L ± ,x a ,x b ,n) 


Q 2 


c l R2 


2 \ -Ff 


P \ V 


= Wfr 12 (L ± ,x a , x b , ^ 

V f 1 

= W?r 12 (-2d fj ,x a ,x b , f x)E?(-2d f i,ri (^)" . 

(4.16) 


We have rewritten the anomaly and the remainder Wj b ^ 12 as functions of the logarithm 
Lj_ = 2 In and have introduced the exponent 

fj=^(C a + C b )T 0 \n%, (4.17) 

47T fl z 

as well as the quantity Ej_(Lj_, fi) which contains the higher-log contributions to the 
anomaly 

E ± (L ± ,v) = exp (~L ± Ff(L ± ) - fjf ± (L ± )) , (4.18) 

and involves the function f±(L±) defined through 

Ff(L ± ) = ^(C a + C b ) T 0 [L ± + 2/_l(L_l)] . (4.19) 

o7r 

The function E±(L±, //) depends on fj, but the derivatives only act on the exponent in 
eq. (4.16). After combining eq. (4.16) with the solution for the two jet functions, which 
involve the evolution factor 


Uj(nj,fx) = exp 


—4(Ci + C 2 )S(nj, fx) + 2 A 1J2 ([Xj, fx) + 2A 7Ja (fXj,fx) , 


(4.20) 
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and the quantity rij. = 2CiA lcusp (/Xj, /x), we obtain the resummed cross section for the 
hadron collider case 
da 

dr±d(cos 8)dx a dxb 

Hjj(Q , 9, //) A*) ^ 12 (-20 j? , 9 , n)E ab (—2d ri , /x) 

_ _ \ f ‘2.cr\0 sin 0 t \ ^ p 'yEV 

^jc 1 ±(d v + hiX j ,n j )j C2± {d v + lnX j ,^ j )—\^— -—- j jyyy, (4.21) 


where 




2/iQ sin 0 

2 5 


V ■= V + Vh + Vj 2 - (4.22) 


To resum all logarithms, one combines eq. (4.21) with the RG evolved hard function from 

[18]. 


5 Two-loop anomalous dimensions for N 2 LL resummation 
5.1 Leptonic collisions 


In this section we compare our results for the lepton-collider case with the fixed-order 
expression for the transverse-thrust spectrum, obtained numerically with the fixed-order 
Monte-Carlo program EVENT2 [37]. If we expand the resummed expression in SCET in 
fixed-order perturbation theory, it has to reproduce the t± —> 0 singularities of the full 
fixed-order result. The EVENT2 code allows us to verify this agreement numerically. More 
importantly, the comparison allows us to numerically determine the two-loop anomalous 
dimensions of the jet and soft functions, the only ingredients that are missing in order to 
achieve N 2 LL accuracy. 

The EVENT2 code provides the 0(a s ) and 0(n 2 ) corrections to differential event-shape 
spectra, i.e. the coefficients A(r ±_) and B(t ±) in the expression 


^=^) + S)^ )+ (£) 2b(tJ+ '"’ p- 1 * 

where a$ is the Born-level cross section. To obtain the transverse-thrust spectrum with 
EVENT2 we modified the usual thrust computation in EVENT2. This can be easily done by 
evaluating the thrust after dropping the longitudinal components of the momenta. 

To obtain the fixed-order expansion of the SCET result, one can set the jet scale fXj, 
soft scale /x s , and hard scale all equal to Q in eq. (4.13). All RG evolution factors 
become trivial in this limit and one obtains 


1 da 
do dr±d(cos 9) 


= lim - (1 + cos 0 2 )H(Q 2 ,Q 2 ) 

. _ 2 p — 'YeV 

jc± {d v , Q) S± (d v , Q ) Tl 1+ ''jyyy ( 4sin2 °y 


(5.2) 


The limit rj —» 0 can be taken after taking the derivatives with respect to it and expanding 
T^_ 1+n in terms of distributions. Using the general expressions for the hard function and 
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Figure 8. Left panel: Comparison of the singular terms Da in eq. (5.3) (solid line) with the one- 
loop coefficient A from EVENT2 (points). Right panel: One-loop term from EVENT2 after subtracting 
the singular terms given by Da- For comparison we show the results both for transverse thrust 
(black) and ordinary thrust (red). 

the Laplace-transformed soft and jet functions, as given, for instance, in ref. [32], with the 
one-loop coefficients and anomalous dimensions from the previous sections, we obtain the 
singular terms in the transverse-thrust distribution. We collect these singular terms in the 
D coefficients, according to 

S(T±)+ (£) Da{t±) + (£) 2 Db(t±)+ ■ ■ ■ ■ (5 - 3) 

where we integrated over 6, to directly compare with the result from EVENT2. The D 
coefficients should reproduce the singularities in t± of the fixed order result in eq. (5.1). 
These coefficients contain plus distributions and 5-functions, but away from rj_ = 0 they 
reduce to regular functions. For t_l / 0, we find 

D a (t l ) = - (-17 + 48 In 2+ 12 In r ± ), (5.4) 

and when we compare with the coefficient A from EVENT2 we find good agreement as t± —>• 
0, as shown in the left panel of figure 8 . To see more clearly that the singular contributions 
are indeed reproduced by the effective theory, we subtract the singular contribution in 
eq. (5.4) from the EVENT2 result, and plot the remainder using logarithmic binning on the 
x-axis. The result of this subtraction is shown in the right panel of figure 8 as black points. 
We can clearly see that the remainder goes to 0 when rj_ —> 0, as it should. For comparison 
the corresponding result for ordinary thrust is also shown in the figure (red points). 

The Db coefficient can be written as 

D b = ^ (M? + D { ^ In t_l + D In 2 r ± + D f In 3 tj_) , (5.5) 

and the one-loop computations in section 3, together with the RG equations for the ingredi¬ 
ents of the factorization formula, determine all coefficients, except for Dg\ which depends 
on the two-loop anomalous dimensions of the soft and jet functions. Because of the relation 
7 // = 7 s j_ + 27 ]^, which stems from RG invariance of the cross section, and because the 
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Cf term 



Figure 9. C\ term in B from EVENT2 after subtracting the Dp\ i = 1,2, 3, singular terms defined 
in eq. (5.5). The horizontal lines corresponds to the C\ part of D^\ The results for transverse 
thrust are shown in black, the ones for ordinary thrust in red color. 

two-loop hard anomalous dimension is known, there is only one unknown coefficient, which 
we take as 7 i s j_. Therefore, we can subtract the singular terms that contain ln*rj_, with 
i = 1, 2,3 from the two-loop results from EVENT2, and the remainder, when multiplied by 
tj_, should be constant when tj_ — > 0. By fitting this constant remainder to the expression 
for D^p in terms of 7 i s _l, we can determine the two-loop soft anomalous dimension. To do 
so, it is useful to separate the three color structures, Cp, CpTpUf, and CfCa , which arise 
at two-loop order. This separation provides an additional check because the part of 7 i s j_ 
proportional to Cp must vanish by the non-abelian exponentiation theorem, since the soft 
function is a matrix element of soft Wilson lines. We show the remainder for the two-loop 
Cp term in figure 9. From the figure, we can see that it is indeed constant, and nicely 
agrees with the Cp part of Dp . For the other color structures, we fit the corresponding 
Dp' 1 terms to obtain 7 i s j_. The results from these fits can be seen in figures 10 and 11, 
where the solid line is the result of a fit to the EVENT2 result, and the band is chosen to 
cover the range spanned by the error bars. The range of lnTj_ where the band and the 
line are plotted corresponds to the fit range. We can see that the remainder in the Cf 
and Ca cases nicely goes to a constant for rj_ —)• 0. There seems to be some numerical 
instability below lnrj_ ~ —11 in the Tp case, and the remainder grows. It is unclear why 
this happens 4 , but the problem also arises for the Tp structure in ordinary thrust, see the 
red points in figure 11. In view of this numerical problem, we perform the fit only for 

4 To obtain our results, we have run EVENT2 in quadruple precision with a value of the CUTOFF parameter 
equal to 10 -18 and the NPOW parameters set to 6. We have also performed runs (up to lnrx = —12) with 
CUT0FF= 10 —12 ,1(T 15 ,10 -16 , and decreasing NPOW to 4, and found similar behaviour. 
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Ca term 



In t_l 

Figure 10. CfCa term in B from EVENT2 after subtracting the Dg\ i = 1,2,3, singular terms 
given in eq. (5.5). The solid line and the band are the result of the fit for the two-loop anomalous 
dimension for transverse thrust (see text). The results for ordinary thrust are shown as the red 
points and dashed line, which corresponds to the known CfCa part of D^\ 
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T F term 



In Tj_ 


Figure 11. CpTp term in B from EVENT2 after subtracting the D^\ i = 1,2,3, singular terms 
given in eq. (5.5) (black points). The solid line and the band are the result of the fit for the two-loop 
anomalous dimension (see text). For comparison results for ordinary thrust are shown as red points, 
together with a dashed line, which corresponds to the known Tp part of D^f 1 for ordinary thrust. 


lnrj_ > —11 in this case. We obtain 

7ls_l_lcu = 148^20 i 71 s_lIt f = 181*12 ’ (5-6) 
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where 


7is-L = : CfCa 7is±| Cj4 + CfTfuj 7i s _l| Tf , (5.7) 

and we recall that, as discussed in previous sections, the above numbers correspond to the 
following definition of the anomalous dimension 

4 kQ sin 2 9 \ f 4 nQ sin 2 9 \ _ / 4 kQ sin 2 9 

In--- ,fj, I = -2 I 7 CUSPCS m-b 7 s ± I I hr- ,/j, 

** ^ ^ (5-8) 

With the value of the two-loop soft anomalous dimension obtained in the equations 
above, and the two-loop jet anomalous dimension determined through the relation 7 ^ = 
( 71 # — 7 is_l)/ 2 , we have now determined all ingredients for N 2 LL resummation accuracy 
for the transverse-thrust differential cross section in leptonic collisions. 

5.2 Hadronic collisions 

To perform N 2 LL resummation in the hadron-collider case, we need the two-loop anoma¬ 
lous dimensions of all the ingredients, together with the two-loop anomaly exponent. The 
general result for the two-loop anomalous dimension of the hard function is known due 
to factorization constraints, and the two-loop jet function anomalous dimension is known 
numerically, from our determination in the previous subsection. This leaves the two-loop 
anomalous dimensions of the beam and soft functions and the anomaly exponent as un¬ 
knowns. However, the soft and beam functions are evaluated at the same scale, so we 
only need the anomalous dimension of their product, i.e. the anomalous dimension of the 
remainder function Wjj in eq. (2.54), which is 

lY™ = 7s_l + 2 7 ? =7 H~ , (5-9) 

see eq. (4.6). The anomaly exponent F± is thus the single unknown ingredient for N 2 LL 
accuracy in the quark-jet case. For gluon jets, also the two-loop gluon-jet anomalous 
dimension 77 * will be needed. We can obtain the two-loop result for 7 ^® by considering 
transverse thrust for the process H —> gg + X, which involves gluon jet functions instead 
of the quark jet functions relevant for 7 * —> qq + X. The two-loop soft functions of the two 
processes are related by Casimir scaling, i.e. the one for the gluon case can be obtained 
by replacing Cf —>• Ca in the result for the qq final state. Since the hard anomalous 
dimensions are also known, 7 ^® follows from RG invariance. Due to Casimir scaling, it is 
also sufficient to determine F± in an arbitrary channel and because the divergences in the 
analytic regulator cancel, a computation of the two-loop 1/a divergence of either the beam 
or the soft function will be sufficient. We will now discuss the simplest way to extract the 
anomaly F±. 

Since soft matrix elements are simpler than collinear ones, it seems preferable to extract 
the anomaly from the soft function. However, it involves Wilson lines along four directions, 
which leads to nontrivial color structure and nontrivial dependence on the scattering angle 
9. Both complications can be avoided by considering transverse thrust for pp -A- Z+ jet and 
pp — > Z + Z. The first one, involving a single Z-boson, is of interest phenomenologically 
and has been measured by CMS [46] . The factorization theorem for this case has the same 
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structure as the one for the two-jet case, but obviously the hard functions are the ones 
relevant for pp —> Z + jet, which were determined to two-loop accuracy in [47]. The jet 
function relevant for the Z-boson is trivial. While there can be hadronic radiation collinear 
to the Z-boson, this effect is power suppressed in t± and does not arise in our leading-order 
factorization theorem. In addition, the soft function involves only three-legs because the 
Z-boson does not carry color charge. The color structure of the soft function is then trivial, 
see e.g. [48]. 

The soft function becomes even simpler in the pp — > Z + Z case, which does not 
involve hard radiation in the final state, except along the beam directions. We can simplify 
things even a bit more by considering energetic electrons instead of Z-bosons, i.e. by 
computing transverse thrust in pp — > 7* — > e + e~ , which involves the same hard function 
we encountered in the lepton-collider case. Neglecting electromagnetic interactions, both 
jet functions are trivial and the only hadronic contributions arise from the beam and soft 
functions. In this case, the soft function will be scaleless and vanishes with the standard 
analytic regulator [22], but we can use the form eq. (2.57) to work with a non-vanishing soft 
function. A simple form of the two-loop soft matrix element can be found in Appendix C 
of [34] . Instead of performing the two-loop computation analytically, one can try to extract 
the anomaly coefficient numerically using a N 2 LO fixed order code for pp —>• 7 * —>• e + e _ , 
such as FEWZ [49] or DYNNLO [50]. To increase the accuracy, one can run this code 
using simple model PDFs, since the anomaly does not depend on their form. 

It is interesting to ask whether this two-loop computation could be avoided by defining 
an e + e _ —>• 2 jets event shape involving the same soft function as transverse thrust in 
pp —> e + e~. For such an observable, one could run EVENT2 and extract the anomaly 
coefficient using the same technique as we used for the jet-function anomalous dimension 
7 ^ in the previous subsection. Naively, it appears that one can achieve this by computing 
transverse thrust along the electron direction in e + e _ —> 2 jets, where transverse is defined 
with respect to the thrust axis. The outgoing jets would play the role of the beams in the 
hadronic collisions and the incoming electron would define the thrust vector. However, an 
interesting complication arises: in contrast to the beam axis, the thrust axis is recoiling 
against soft radiation. This changes the factorization theorem, which then has the same 
form as the one for jet broadening [23, 24]. In the presence of recoil, the soft function is 
not exactly the same as the one shown in eq. ( 2 . 21 ), but will have to be computed at a 
fixed value of the transverse momentum in each hemisphere (where the hemispheres are 
defined by the thrust vector). The soft transverse momentum is opposite and equal to 
the transverse momentum of the collinear radiation in the given hemisphere and the event 
shape will then involve an integral over the transverse momentum. To avoid this problem 
and make the e + e _ —>• 2 jets case more similar to the pp —> e + e~ case, one could try 
to use a recoil-free definition of the jet axis, such as the broadening axis or the winner- 
take-all axis [51]. This solves the recoil problem, but the two axes will in general not be 
back to back, and in order to define them, one will first need to split the event into two 
hemispheres, for example using the thrust axis. The resulting soft function is thus again 
more complicated than the original one, and the two-loop anomaly coefficient will likely 
be a function of the angle between the thrust axis and the broadening axis and extracting 
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this function numerically seems difficult. We conclude that a numerical extraction using 
pp 7* ~► e+ e_ looks more promising. We will not pursue the extraction further in the 
present paper, but plan to come back to this issue in the future when we implement the 
resummed expression for transverse thrust numerically. 

6 Conclusions and outlook 

In this paper we have analyzed transverse thrust in the dijet limit. Our findings are 
synthesized in eq. (2.55), the factorized expression of the transverse-thrust cross section for 
hadronic colliders in Laplace space. This result, which we derived within the framework of 
SCET, provides the basis for all-order resummations of enhanced perturbative corrections 
to this observable, beyond the NLL accuracy which has been achieved in the literature. 

The factorization formula for transverse thrust is quite interesting from the point of 
view of the effective theory and displays several remarkable features, which are worth em¬ 
phasizing. In particular, it involves collinear modes at more than one invariant mass scale, 
together with soft modes, and therefore contains the two regimes of SCET, called SCETi 
and SCETn, together in the same problem; this is the first collider-physics example that 
we are aware of, where this occurs. It also involves a collinear factorization anomaly, beam 
and jet functions, as well as matrix-valued soft and hard functions, bringing together in one 
problem several different effective-theory objects that were developed and studied in re¬ 
cent years. Anomaly cancellation involves an intricate interplay between the soft and beam 
functions, which leads to all-order constraints on the form of the collinear anomaly, similar 
to the constraints factorization imposes on the structure of soft anomalous dimensions. 

As an instructive starting point of our studies, we analyzed transverse thrust in leptonic 
collisions, which provide a simplified environment to study transverse event shapes. The 
resulting factorization formula is much simpler than the hadron-collider result and has the 
same structure as the one for ordinary thrust. Using our result, we were able to numerically 
extract a two-loop anomalous dimension that is needed in the hadron-collider case from 
an analysis of the leptonic result; this provides a nice example of the universality which 
becomes manifest after separating physics associated with the different relevant energy 
scales. 

In order to obtain N 2 LL resummation accuracy for hadronic colliders from our general 
factorization formula, one needs the two-loop anomalous dimensions as well as the two- 
loop anomaly exponent. Using RG invariance and the general results for the anomalous 
dimensions of the hard function, together with our result for the lepton-collider case, we 
obtain all two-loop anomalous dimensions. The single missing ingredient for N 2 LL accuracy 
is thus the two-loop anomaly coefficient, which can be extracted numerically by computing 
transverse thrust in Drell-Yan production using one of the existing N 2 LO fixed-order codes. 

In order to perform a phenomenological analysis and compare with data, one needs 
to match the resummed results to a fixed-order computation [45], and to carefully study 
power corrections. Our results for leptonic collisions seem to indicate that the matching 
to fixed-order may be a larger correction than for regular thrust. In the future, we will 
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extract the missing ingredient for hadronic collisions numerically and implement the N 2 LL 
resummed expression for transverse thrust into a numerical code. 

It would be interesting to extend the present analysis to other hadronic event shapes. 
Instead of analytically computing the ingredients of the associated factorization theorems 
for each given observable, it would be much more efficient to perform a fully numeri¬ 
cal evaluation of the corresponding soft, jet, and beam functions; this would provide an 
effective-theory-based, automated framework for computing such observables, which could 
be integrated with efforts to construct a SCET-based Monte Carlo event generator [53]. 
There exists already an automated resummation framework, CAESAR [8], for hadronic 
event shapes at NLL accuracy, and for a class of lepton-collider event shapes an exten¬ 
sion to N 2 LL was recently achieved in ref. [54]. Also within SCET, automated N 2 LL 
resummation has been performed, so far not for event shapes, but for cross sections with 
electroweak bosons in the presence of a jet veto [55], using the automated fixed-order NLO 
code MadGraph5_aMC@NLO [56]. 

An interesting application of the event-shape computations is that they could provide 
an alternative subtraction scheme for dijet and H/Z + jet production, along the line of the 
qx subtraction scheme by Catani and Grazzini [57]. The resummed cross section includes 
the virtual corrections as well as the singularities arising for t± —> 0. It can thus be used as 
a subtraction in this limit. Configurations with t_l > 0, on the other hand can be computed 
using the NLO prediction for Z + 2jets (or 3jets for dijet observables). 

Before closing, let us mention that we did not consider Glauber gluons in the derivation 
of our factorization formula. Our result assumes the standard factorization of the cross 
section into a perturbative kernel convolved with PDFs and will resum the logarithms en¬ 
countered in the perturbative computation of the hard-scattering kernel. Recently ref. [58] 
showed that the standard diagrammatic arguments (see e.g. refs. [59-62]) to show the 
absence of Glauber contributions fail for observables such as transverse thrust. It would 
be interesting to analyze the contribution from Glauber gluons within SCET but to date 
a complete implementation of this mode into the effective theory is not yet available. If 
Glauber modes are indeed present in transverse thrust, they could mediate spectator- 
interaction contributions to the cross section, in addition to the terms captured by our 
result. These would likely involve non-perturbative physics not encoded in the PDFs. A 
comparison of our factorized results with data may help shed some light on these effects 
and how they relate to the underlying-event contribution supplied by parton-shower Monte 
Carlo programs. Once their form is understood, our result can also be used to study these 
effects quantitatively. To do so, one could use a combination of event shapes that maxi¬ 
mizes the sensitivity to underlying-event effects. On the other hand, to, for instance, obtain 
a determination of the strong coupling from a hadronic event shape analysis, one would 
choose a combination that is as less sensitive as possible to underlying-event effects. We 
look forward to investigate these issues in the future. 
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A Anomalous dimensions 


The coefficients of the cusp anomalous dimension that are needed for N 2 LL accuracy are 
given by 


r n = 4. 
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is the QCD beta function. The quantity T 7j {v, /i) is defined like T 7 cusp (zv, fi) but with 7 CU sp 
replaced by 7 j i . 

The coefficients of the quark and gluon anomalous dimensions y 9 and 'y 9 , respectively, 

read 
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The color factors are given by 

Cf=N wT’ Ca = N c , T f = \, (A.11) 

where N c is the number of colors, and nj the number of light fermions. 


B One-loop soft, jet, and beam functions 
B.l Jet functions 


In the main text, we gave the result for the quark jet function in light-cone gauge. For 
convenience we now also give the result of the diagrams in Feynman gauge, where the first 
three diagrams in figure 5 contribute. The contribution from the first diagram is 


Dij(r c ±) = a s C F 


(27t) 3 2e r c _l V 

the second and third one are identical and give 
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where = 2 F d ^ 2 /T(d/2) is the d-dimensional solid angle. The reader can easily check 
that the sum of the diagrams gives the same result which is obtained in light-cone gauge 
in eq. (3.3). 

In light-cone gauge, the bare gluon jet function is obtained by computing the fifth and 
sixth diagrams in figure 5 


Jg± e ( T c±) = S(tc±) + D 5j + D ej , 


for which we obtain 
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B.2 Beam functions 


(B.3) 


(B.4) 

(B.5) 


The diagrams that contribute to the beam functions at one loop are shown in figure 6. We 
denote the contribution of the zth diagram in the figure by Di F when it corresponds to the 
c a sector, and by Die when it corresponds to the Cb sector. 

The diagrams in the first row of figure 6 contribute to the quark beam function. In 
Feynman gauge, the first diagram does not need the analytic regulator to be well defined 
and, therefore, gives the same result in the c a and Cb sectors. We obtain 

Di B (z,t ± ) = D 1b (z,t±) = ^^ £ C'(e)(l -e)(l - z) — (— U —~\ , (B.6) 
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with V'l 31 ) the digamma function, and it is understood that z corresponds to x a (x b ) in 
the c a (cfe) sector. For the sum of the second and third diagrams, which give identical 
contributions, we obtain 
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where <f>(, 2 , s, a) is the Lerch transcendent. The fourth diagram in the first row vanishes. 

The diagrams in the second row of figure 6 contribute to the gluon beam function, and 
its computation is analogous to the quark case above. We obtain 
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and the eighth diagram vanishes. 

The diagrams in the third row of figure 6 contribute to the off-diagonal coefficients 
Iq-^g, first diagram in the row, and / 9 <_ 9 , second diagram in the row. We obtain 

D^rA = A»(.,rJ = (l - i (^,) 2 '. (B-14) 
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The divergences in the analytic regulator for all the expressions above can be made 

manifest by expanding 
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B.3 Soft functions 


The X integrals in eq. (3.30) are given by 
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where it is understood that we take it at a = 0 for the I 12 integral, which does not require 
the analytic regulator. For convenience, we also define 
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For arbitrary values of a , one has Z' + {6) = X + (—#) and Z'_(Q) = —X_(—#). The numerical 
values of the constants are H ~ —1.85939 and XT ~ 8.44015. The results for the func¬ 
tions a(B(#) and #B)(#) are plotted in figure 12 , as blue solid and magenta dashed lines, 
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Figure 12. Numerical evaluation of the a^(6) (blue solid line) and &W(0) (magenta dashed line) 
terms in eqs. (B.26) and (B.27). 


respectively. The functions are symmetric in 9. For positive 6 , they are obtained from a 
numerical evaluation of the expressions 


a (1) (0) = 


f dx ( ,5( x ) (i n 2 -g 2 In sin 9) 2 

I o 8 tt(2 - x) \ 


+ 


x Im 


(4 In 2 — 2 ln(2 — x) + 4 In sin 9) + 

1 , (yi-V 1 + yi\\ 


lnx 


J + 


{ V 1 + vi V^ 1 + V 1 + y‘i, 


(B.29) 


6(1) (0) = L dx ( 6{x) ( 3 ln 2+2 ln cot 0 


+ 


l*J + 


Q 

8 In 2 — 2 ln(2 — x) + 4 ln cot - ) + 


lnx 


f (1 - cos6») (1 + yi(l - x)sm9 + yj(l + cos6»)) f Vl 
x 1m < --- , -- ln — 


-Vi + vi 


\/i + y( 


+ V± + vI 


y i 


(B.30) 


where y\ was given in eq. (3.36). 
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C Lepton tensor 

The lepton tensor, including contributions from photon and Z exchanges, is given by 


T V = _ 

^ 2Q 2 


= : -L v ( g, w - 2 


9fj. v 2 


V\txV2v + P2,xP\v 

Q 2 


2 2Q 2 v e v f Qf Q 4 (v 2 + aj)vj 
Q 2 ~M 2 Z + ( Q 2 -M 2 ) 2 


^Pl/j,P2u + P2nPlv 

Q 2 


(C.l) 


T^= — 

IMJ 2Q 2 


4 7 _ 0 Pljj,P2v + P2nPlv\ Q A ( v l + a l) a ) 

9i-iu 2 


=■ ~L [g^-2 


Q 2 J (Q 2 - M|)2 

PlnP2v +P2^,Plu' 

Q 2 


(C.2) 


with pi and p 2 the electron and positron momenta, respectively. The vector and axial 
charges, Vf and a/, are defined as 


v f : = 


Tj — 2QfS 


2 swcw 


rp 3 
W 

; a f := 


2.S]yCw ’ 


(C.3) 


with Tj the third component of the weak isospin, Qf the electric charge of fermion / (with 
Qe = 2T e 3 = -1), cw : = cos 9wi sw ■= sin^vT) the weak mixing angle, and Mz the 
Z-boson mass. 


References 

[1] A. Banff G. P. Salam and G. Zanderighi, JHEP 0408 (2004) 062 [hep-ph/0407287], 

[2] A. Banff, G. P. Salam and G. Zanderighi, JHEP 1006 (2010) 038 [arXiv: 1001.4082 [hep-ph]]. 

[3] V. Khachatryan et al. [CMS Collaboration], Phys. Lett. B 699 (2011) 48 [arXiv:1102.0068 
[hep-ex]]. 

[4] G. A ad et al. [ATLAS Collaboration], Eur. Phys. J. C 72 (2012) 2211 [arXiv:1206.2135 
[hep-ex]]. 

[5] G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 88 (2013) 3, 032004 [arXiv:1207.6915 
[hep-ex]]. 

[ 6 ] V. Khachatryan et al. [CMS Collaboration], JHEP 1410 (2014) 87 [arXiv: 1407.2856 [hep-ex]]. 

[7] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. D 83 (2011) 112007 [arXiv:1103.5143 
[hep-ex]]. 

[ 8 ] A. Banff, G. P. Salam and G. Zanderighi, JHEP 0503 (2005) 073 [hep-ph/0407286], 

[9] C. W. Bauer, S. Fleming, D. Pirjol and I. W. Stewart, Phys. Rev. D 63 (2001) 114020 
[hep-ph/0011336]. 

[10] C. W. Bauer, D. Pirjol and I. W. Stewart, Phys. Rev. D 65 (2002) 054022 [hep-ph/0109045], 

[11] M. Beneke, A. P. Chapovsky, M. Diehl and T. Feldmann, Nucl. Phys. B 643 (2002) 431 
[hep-ph/0206152], 

[12] T. Becher, A. Broggio and A. Ferroglia, arXiv: 1410.1892 [hep-ph]. 


- 45 - 



[13] T. Becher and M. Neubert, Eur. Phys. J. C 71 (2011) 1665 [arXiv: 1007.4005 [hep-ph]]. 

[14] J. C. Collins, D. E. Soper and G. F. Sterman, Nucl. Pliys. B 250 (1985) 199. 

[15] V. Mateu and G. Rodrigo, JHEP 1311 (2013) 030 [arXiv:1307.3513], 

[16] S. Etter, MSc thesis. University of Bern (2014). 

[17] I. W. Stewart, F. J. Tackmann and W. J. Waalewijn, Phys. Rev. D 81 (2010) 094035 
[arXiv:0910.0467 [hep-ph]]. 

[18] R, Kelley and M. D. Schwartz, Phys. Rev. D 83 (2011) 045022 [arXiv: 1008.2759 [hep-ph]]. 

[19] A. Broggio, A. Ferroglia, B. D. Pecjak and Z. Zhang, JHEP 1412 (2014) 005 
[arXiv:1409.5294 [hep-ph]]. 

[20] T. Becher and M. Neubert, JHEP 0906, 081 (2009) [Erratum-ibid. 1311, 024 (2013)] 
[arXiv:0903.1126 [hep-ph]]. 

[21] J. -y. Chiu, F. Golf, R. Kelley and A. V. Manohar, Phys. Rev. D 77 (2008) 053004 
[arXiv:0712.0396 [hep-ph]]. 

[22] T. Becher and G. Bell, Phys. Lett. B 713 (2012) 41 [arXiv:1112.3907 [hep-ph]]. 

[23] J. -y. Chiu, A. Jain, D. Neill and I. Z. Rotlrstein, Phys. Rev. Lett. 108 (2012) 151601 
[arXiv:1104.0881 [hep-ph]]. 

[24] T. Becher, G. Bell and M. Neubert, Phys. Lett. B 704 (2011) 276 [arXiv: 1104.4108 [hep-ph]]. 

[25] T. Becher and M. Neubert, JHEP 1207 (2012) 108 [arXiv:1205.3806 [hep-ph]]. 

[26] T. Becher, M. Neubert and L. Rothen, JHEP 1310, 125 (2013) [arXiv: 1307.0025 [hep-ph]]. 

[27] M. G. Echevarria, A. Idilbi and I. Scimemi, Phys. Lett. B 726 (2013) 795 [arXiv:1211.1947 
[hep-ph]]. 

[28] S. Catani and M. H. Seymour, Phys. Lett. B 378 (1996) 287 [hep-ph/9602277]. 

[29] T. Becher and M. Neubert, Phys. Rev. Lett. 102, 162001 (2009) [Erratum-ibid. Ill, no. 19, 

199905 (2013)] [arXiv:0901.0722 [hep-ph]]. 

[30] E. Gardi and L. Magnea, JHEP 0903, 079 (2009) [arXiv:0901.1091 [hep-ph]]. 

[31] L. J. Dixon, E. Gardi and L. Magnea, JHEP 1002, 081 (2010) [arXiv:0910.3653 [hep-ph]]. 

[32] T. Becher and M. D. Schwartz, JHEP 0807 (2008) 034 [arXiv:0803.0342 [hep-ph]]. 

[33] T. Becher, M. Neubert and B. D. Pecjak, JHEP 0701 (2007) 076 [hep-pli/0607228]. 

[34] T. Becher and G. Bell, JHEP 1211 (2012) 126 [arXiv:1210.0580 [hep-ph]]. 

[35] T. Becher and G. Bell, Phys. Lett. B 695, 252 (2011) [arXiv: 1008.1936 [hep-ph]]. 

[36] J. C. Collins and D. E. Soper, Nucl. Phys. B 194, 445 (1982). 

[37] S. Catani and M. H. Seymour, Nucl. Phys. B 485 (1997) 291 [Erratum-ibid. B 510 (1998) 
503] [hep-ph/9605323]. 

[38] R. Kelley and M. D. Schwartz, Phys. Rev. D 83 (2011) 033001 [arXiv: 1008.4355 [hep-ph]]. 

[39] E. W. N. Glover and M. E. Tejeda-Yeomans, JHEP 0306, 033 (2003) [hep-ph/0304169]. 

[40] E. W. N. Glover, JHEP 0404, 021 (2004) [hep-ph/0401119]. 

[41] Z. Bern, A. De Freitas and L. J. Dixon, JHEP 0306, 028 (2003) [Erratum-ibid. 1404, 112 
(2014)] [hep-ph/0304168]. 


- 46 - 



[42] A. De Freitas and Z. Bern, JHEP 0409, 039 (2004) [hep-ph/0409007, hep-pli/0409007]. 

[43] Z. Bern, A. De Freitas and L. J. Dixon, JHEP 0203, 018 (2002) [hep-ph/0201161]. 

[44] T. Becher and M. Neubert, Phys. Rev. Lett. 97, 082001 (2006) [hep-ph/0605050]. 

[45] Z. Nagy, Phys. Rev. D 68 (2003) 094002 [hep-ph/0307268]. 

[46] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 722, 238 (2013) [arXiv: 1301.1646 
[hep-ex]]. 

[47] T. Becher, G. Bell, C. Lorentzen and S. Marti, JHEP 1402, 004 (2014) [arXiv: 1309.3245 
[hep-ph]]. 

[48] T. Becher, G. Bell and S. Marti, JHEP 1204, 034 (2012) [arXiv:1201.5572 [hep-ph]]. 

[49] R. Gavin, Y. Li, F. Petricllo and S. Quackenbush, Comput. Phys. Commun. 182, 2388 
(2011) [arXiv:1011.3540 [hep-ph]]. 

[50] S. Catani, L. Cieri, G. Ferrera, D. de Florian and M. Grazzini, Phys. Rev. Lett. 103, 082001 
(2009) [arXiv:0903.2120 [hep-ph]]. 

[51] A. J. Larkoski, D. Neill and J. Thaler, JHEP 1404, 017 (2014) [arXiv:1401.2158 [hep-ph]]. 

[52] M. Grazzini, JHEP 0802, 043 (2008) [arXiv:0801.3232 [hep-ph]]. 

[53] S. Alioli, C. W. Bauer, C. J. Berggren, A. Hornig, F. J. Tackmann, C. K. Vermilion, 

J. R. Walsh and S. Zuberi, JHEP 1309 (2013) 120 [arXiv:1211.7049 [hep-ph]]. 

[54] A. Banfi, H. McAslan, P. F. Monni and G. Zanderighi, arXiv:1412.2126 [hep-ph]. 

[55] T. Becher, R. Frederix, M. Neubert and L. Rothen, arXiv:1412.8408 [hep-ph], 

[56] J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer, H.-S. Shao and 
T. Stelzer et al. , JHEP 1407, 079 (2014) [arXiv:1405.0301 [hep-ph]]. 

[57] S. Catani and M. Grazzini, Phys. Rev. Lett. 98, 222002 (2007) [hep-ph/0703012], 

[58] J. R. Gaunt, JHEP 1407 (2014) 110 [arXiv:1405.2080 [hep-ph]]. 

[59] G. T. Bodwin, Phys. Rev. D 31 (1985) 2616 [Erratum-ibid. D 34 (1986) 3932], 

[60] J. C. Collins, D. E. Soper and G. F. Sterman, Nucl. Phys. B 261 (1985) 104. 

[61] J. C. Collins, D. E. Soper and G. F. Sterman, Nucl. Phys. B 308 (1988) 833. 

[62] J. Collins, “Foundations of perturbative QCD,” (Cambridge monographs on particle physics, 
nuclear physics and cosmology. 32) 


- 47 - 



